Ore A AT A 


WOl 2. NO°4 


December 1955 


Reflection by semi-infinite diffusers 


R. G. GIOVANELLI 


Division of Physics, National Standards Laboratory, 
Commonwealth Scientific and Industrial Research 
Organization, Sydney. 


SUMMARY. — Total and directional reflectances are derived and tabulated for semi-infinite diffusers with matrices of various refractive 
indices and scattering according to the phase function o(1 + 2x cos 0). 


SOMMAIRE. — On établit les formules algébriques et on donne des lables des valeurs numériques des réflectances totales et direct ion- 
nelles, calculees pour des difjuseurs semi-infinis caractérisés par des indices et des coefficients de diffusion satisfaisant ala loi 


Gp (1 + x cos 8). 


ZUSAMMENFASSUNG. 


Der Gesamwert und der gerichlete Anteil der Ritckstrahlung eines den unendlichen Halbraum erfiillenden 


streuenden Mediums werden fiir Matrizzen von verschiedenen Brechwerten und einer Streuung gemdss der Phasenfunktion 


Dp (1 + x cos 0) berechnet und tabelliert. 


I. Introduction. — There are potentially at least 
two major applications of a precise knowledge of the 
total and directional reflectances of diffusing mate- 
rials, namely colour prediction and spectro-photome- 
tric chemical analysis. 

A number of attempts have been made to describe 
the reflectance and transmittance of plane parallel 
diffusers in terms of the scattering and absorbing 
properties of their constituents; by KuBELKA and 
Munk [1], Amy [2] and KusBeLxa [3] using the 
SCHWARZCHILD approximation in solving the relevant 
equations of radiative transfer, and by others. Prrrs [4] 
has studied the light scattered in photographic emul- 
sions using the Enpinaton and Mark approximate 
methods. Exact methods have been developed by 
CHANDRASEKHAR [5] for obtaining some of the optical 
properties of diffusers, including the reflected intensity 
distribution for collimated light incident on sem1-infi- 
- nite diffusers scattering isotropically, according to a 
phase function m,(1 + x cos 6) or to the RAYLEIGH 
phase function ; the two former cases include absorp- 
tion within the medium, and tables of some of the 
relevant H-functions needed in evaluating the solu- 
tions have been given for isotropic and RayLeicu 
phase functions, and for the phase function 
wo(1 + cos 9). 

The above mentioned researches have discussed 
primarily diffusers whose matrices have unit refractive 
index. Whenever allowance has been made for matrices 
of higher refractive index, in which case surface 
reflections occur internally, the problem has been sim- 
plified by assuming a reflection coefficient indepen- 
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dent of direction ; 1. e. the surface has been assumed 
perfectly diffusing for internally incident radiation 
reflected and transmitted through it. 

The present paper discusses total and directional 
reflectances of semi-infinite diffusers of constant com- 
position, scattering according to a phase function of 
the form m,(1 + a cos 0). Part I deals with media 
with N (matrix) = 1. Part Il deals with media in 
which NV (matrix) exceeds unity. Throughout, an 
accuracy has been sought better than at present obtain- 
able in best measurements, 1. e. of the order of 
+ 0.001 or better in reflectance. 


Part I. Marrices oF Rerractive InpEx Unity 


II. Isotropic scattering. — Let a collimated beam 
of radiation of flux 7Fi(u»9) per unit area normal to 
the beam be incident on the surface of a semi-infinite 
diffuser at an angle to the outward normal whose 
cosine is —u,» Then CHANDRASEKHAR [5], (§ 26) has 
shown that for isotropic scattering the reflected inten- 
sity in a direction yu. is 


oe ira, 


where zp, the albedo for single scattering, is defined by 
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so and « being the scattering and absorption coeffi- 
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cients respectively, while H(y) is a function defined by 
the equation 


ae) 


(2.3) Hy) =1+ ut) | —*4 Ay! ay’. 
Jo e+e 
; F 1 oe 
For isotropic scattering, W(u’) = ay and the funce- 


tion H(u) has been tabulated by CHANDRASEKHAR [5], 
(Table XI). 
The reflected flux per unit area of surface is 


[S 
iN 
ad 
=| 
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d2 being an element of solid angle, while the total 
reflectance for light incident in a direction yo is given 
by 


(2.5) Ru) = 


It follows readily from the properties of H- functions 
that 


(2.6) Ro) = 1 — H(0) (1 — m9)". 

If the incident radiation is uniformly diffused, with 
flux f dQ per unit area normal to the beam in an ele- 
ment of solid angle dQ, then from (2.1) the reflected 
intensity in a direction yu is 


(2.7) = fA aot 
a result which could also have been obtained directly 
from (2.6) from a well known photometric reciprocity 
law. 

The total reflected flux per unit area of surface is 


fd 
nF =2eff a [—H) (mol) dy 
sh tt) 


74 
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The reflectance for totally diffused incident radiation 
is thus 

m4 

Qo) Re 196 att | bisheneka 
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The value of the integral, the first moment of A(u), 
is tabulated by Cuanpraseknar [5],(Table XX XIII). 


III. Scattering according to a phase function 
mo(1 + x% cos 0). — Scattering by real media rarely 
occurs isotropically. The angular distribution of scat- 
tered light may be described by the phase function 
p(cos 8), where p(cos 0) dQ/4 mm» is the fraction of 
scattered radiation entering solid angle dQ : thus 


[ Orr. Acta 


— Wo ~ 


fe da 
(3.1) | p(cos 9) ace 


e 


For a phase function m,(1 + x cos 0), the reflected 
intensity J(u, e) of a semi-infinite medium in direc- 
tion vw, azimuth 9, is given by CHANDRASEKHAR [9], 
(§ 46) 


(3.2) 


| ‘ 
I(u, ¢) = We Bo Filo) X 


o) U. ol 


 H(u) Ho) [1 e(w + 40) 21 


+ 2(4—p2)2(1 — 2)? Hy) HM (U9) cos(o— @) 
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the direction and azimuth angle of the incident colh- 
mated beam being — yy and 9, respectively. H(y) and 
Hu) are H-functions defined by (2.3), with 


and 
1 
¥(u) = ae, t9(1 — y?) 


respectively, while 
2S 


c = x(1 — wo) mo a 


Df 
™o %o 
a and «, being the moments of H(u) of order zero and 
unity : 
74 74 
to = | A(y) dy, op = | Ay) vu dy. 
0 0 


Integration around ¢ yields for the element of re- 
flected flux in du per unit area of surface 


FW Li(9) U Uo 
2(u + po) 


[£— elu + Yo) —2(1— mo) vo] duo - 


7d, = H(u) H(u9) X 


Hence the total reflected flux is 
(3 a) ni, = 
™o 


= mHi(u9) Uo E — [T(u9) 1— a (% — ¢ a) . 
and the total reflectance for light incident in a direc- 
tion po 18 

(3.4) 


R (uo) = fl (aig) 


wo \ 1 
een (% — € a) \ . 
If the incident radiation is uniformly diffused, the 
reflected flux per unit area is 


a 


(3.5) aFy= | R(u) vo f a2 = 


= nf[1—20,)4—5 (%—ea){ |. 
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Thus the reflectance for totally diffused incident 
radiation is 


(3.6) Ro = 1— 2a, (heat OG ec iay 


Values of the relevant H-functions, of the moments 
% and «, and of the constant c are given by CHAN- 
DRASEKHAR [5], (Tables XVI and XVII) for x = 1. 


IV. Scattering according to a phase function 
mol + a cos 0) ; Eddington’s approximation. — 
The exact results given in § III lead to numerical 
values for the reflectance Eyes tables of relevant 
H-functions are available ; these are only at present 
available for the cases x = 1 and x = 0 (§ I). Using 
EDDINGTON’s approximation, however, reflectances 
may readily be obtained for any value of x. While 
these differ shghtly from the exact reflectances for 
x =O and x = 1, they can be used to provide very 
good interpolations for other values of «. 

Pirts [4] has discussed the behaviour of plane pa- 
rallel diffusing media for light incident normally on 
the surface, and has shown that if the phase function 
be expanded in a series of LEGENDRE polynomials, the 
total intensity is everywhere independent of terms of 
the second and higher degree, on the EppINGToN 
approximation. Thus to this approximation the phase 
function is fully described by m,(1 + 2 cos 6). 

Pirrs’s discussion can be extended in a straight- 
forward manner to the case of radiation incident on a 
semi-infinite diffuser in a direction py, and it can be 
shown that the total reflectance is 


' oat 
(4.4) R(g,) =——= op Tne 
2Vx+38—m 2 | L+povz J 
where 


7, = (3 — mm 2)(1 — m). 


The reflectance for totally diffused incident radia- 
tion follows readily by integration 


Mo \ 


a+ 


(4.2) i= = 
IN) fe 


} 6 a 2(1 — wo) a [Va y — log(1 +  \/ x) ; 
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In the case of isotropic scattering, # = 0 and (4.1) 


and (4.2) simplify to 


(4.3) Ru) = = as 
(142 ve (1+ 32) 
and 
easy: = Peel S2— ats +V39)] 
2 3 4) 


where } = 1 — mp. 
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TABLE 1 


Total Reflectances for isotropic scattering and for a phase 
function @, (1 + cos 6) 


ee. Diffused incident 
= 0 aa Radiation 
%0 
Isotropic | @p(1+ cos 8) Isotropic | G(1-+ cos 8) 

1.000 1.00000 1.00000 1.00000 1.00000 
0.999 0.91285 0.89367 0.92971 0.91446 
0.995 0.81705 0.77877 0.84985 0.81945 
0.990 0.75275 0.70270 0.79457 0.75482 
0.975 0.64092 0.57344 0.69501 0.64140 
0.950 0.53555 0.45552 0.59667 0.53311 
0.925 0.46655 0.38104 0.52965 0.46172 
0.900 0.41495 Ops 2 72 0.47802 0.40825 
0.85 0.33966 — 0.40017 — 
0.80 0.28526 0.20015 0.34187 0.27406 
0.7 0.20867 0.13286 0.25655 0.19626 
0.6 0.15541 0.09065 0.19471 0.14318 
0.5 0.11521 0.06192 0.14653 0.10411 
0.4 0.08336 0.04147 0.10934 0.07394 
0.3 0.05721 0.02638 0.07445 0.04986 
0.2 0.03524 0.01513 0.04626 0.03018 
0.1 0.01639 0.00649 0.02170 0.01382 
0 0.00000 0.00000 0.00000 0.00000 
V. Results. — Table 1 contains data, computed 


from the exact results of sections II and III, on the 
total reflectances for radiation incident along the nor- 
mal or perfectly diffused, for the cases of isotropic 
scattering and scattering according to the phase func- 
tion m»(1 + cos 6). esate are given for values of wo 
tabulated by CHANDRASEKHAR, and for some addi- 
tional values of mp, close to unity for which supple- 
mentary H-functions (not obtainable to sufficient 
accuracy by interpolation) have been computed. 


TABLE 2 


Reflectances obtained with the EppINGTON approximation 
for isotropic scattering and fora phase function a,(1 + cos 6) 


Uy —— al 
Isotropic scattering @o(1 + cos 0) 
Mp 

R(1) Correction R(1) Correction 
1.000 1.00000 0.00000 1.00000 0.00000 
0.999 0.91376 —().00091 0.89484 —0.00117 
0.995 0.81952 —0.00247 0.78125 —0(0.00247 
0.990 0.75649 —().00374 0.70628 —0.00358 
0.975 0.64722 —0.00630 0.57885 —0.00541 
0.950 0.54426 —().00871 0.46270 —0.00718 
0.925 0.47666 —0(0.01011 0.38922 —0.00818 
0.900 0.42596 —0.01101 0.33590 —0.00878 
0.85 0.35152 —0.01186 0.26104 _ 
0.80 0.29729 —0.01203 0.20957 —0.00942 
On 0.22005 —().01138 0.14183 —0.00897 
0.6 0.16548 —0(0.01007 0.09873 —().00808 
0.5 0.12372 —0.00851 0.06888 —0.00696 
0.4 0.09017 —0(0.00681 0.04713 —0.00566 
0.3 0.06230 —().00511 0.03072 —0.00434 
0.2 0.03860 —0(0.00336 0.01802 —0.00289 
0.1 0.01806 —0(.00167 0.00801 —0.00152 
0 0.00000 0.00000 0.00000 0.00000 


Table 2 lists total reflectances, computed for nor- 
mal incidence from the approximate equation (4.1), 
for « = 0 (isotropic scattering) and x = 1. The correc- 
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tions to the two latter sets of results may be derived 
by comparison with columns 2 and 3 of Table 1, and 
are listed in Table 2, from which it will be seen that 
the two sets of corrections are close together, diffe- 
ring at most by little more than 0.002. Consequently 
the corrections to results derived from (4.1) for inter- 
mediate values of x can be obtained reliably by linear 


TABLE 3 
Reflectance for a phase function a, (1 + x cos 0) 
oy == 
x 
eX) 

0 0.25 0.5 0.75 1.00 
1.000 1.00000 1.00000 1.00000 1.00000 1.00000 
0.999 | 0.91285 | 0.9090 0.9046 0.8996 | 0.89367 
0.995 0.81705 0.8094 0.8006 0.7906 0.77877 
0.990 | 0.75275 | 0.7427 0.7312 6.7180 0.70270 
0.975 | 0.64092 | 0.6273 0.6116 | 0.5940 0.57344 
0.950 | 0.53555 | 0.5191 0.5006 0.4796 0.45552 
0.925 | 0.46655 | 0.4488 0.4289 0.4075 0.38104 
0.900 | 0.41495 | 0.3965 0.3760 0.3531 0.32712 
0.85 0.33966 | 0.3209 0.3003 0.2772 0.2516 
0.80 0.28526 | 0.2668 0.2468 0.2246 0.20015 
0.7 0.20867 | 0.1919 | 0.1738 0.1542 0.13286 
0.6 0.15541 | 0.1408 Onmi252 0.1085 0.09094 
0.5 0.11521 | 0.1029 0.0900 0.0764 0.06192 
0.4 0.08336 | 0.0735 0.0633 0.0526 0.04147 
0.3 0.05721 | 0.0498 0.0423 0.0344 0.02638 
(1 | 0.03524 | 0.0304 0.0254 0.0203 0.01513 
0.1 0.01639 | 0.0140 0.0115 0.0106 0.00649 
(0) 0.00000 0.00000 0.00000 0.00000 0.00000 


interpolation. In Table 3 total reflectances for nor- 
mally incident radiation obtained in this way are 
listed for x = 0, 0.25, 0.50, 0.75 and 1.00; reflec- 
tances for intermediate values of 2 may be obtained 
correct to one unit in the third decimal place by linear 
interpolation. The reliability of the above procedure 
has been checked by evaluating a few H-functions 
for x = 0.5 and comparing the reflectances so derived 
with the corrected values. 


Part. I]. Marrices of REFRAcTIVE INDEX 
ExcreEepInG UNItTy 


VI. Introduction. — The refractive index of the 
matrix of a diffuser affects its reflecting properties 
because not all the radiation incident on the surface 
from within passes through, but some is internally 
reflected. Hence for a given albedo for single scatte- 
ring and a given phase function the reflectance still 
depends on the refractive index and on the state of 
the surface. The angular distribution of light after 
internal reflection at the surface has, of course, an 
infinite range of possibilities when all types of surface 
are considered, but fortunately in many cases of inte- 
rest the surface may be regarded as a specular reflec- 
tor, a case discussed here in some detail. An easier 
case to discuss is that of a hypothetical surface which 
diffuses perfectly all radiation reflected from or trans- 
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mitted through it ; its importance is that it provides a 
limiting case, and that the total reflectances of all real 
media must lie in the range established by these two 
cases. 

Results for media with specularly reflecting sur- 
faces and phase functions of the form mo(1 + x cos 8) 
depend on the solution of integral equations, and 
the labour involved becomes exceedingly great 
as the desired degree of approximation increases. An 
approximate formula can be derived, however, for 
the reflectance, using EppiINGTon’s approximation, 
and this may be used to interpolate with adequate 
accuracy between various values of x, thus reducing 
considerably the effort in securing results. 

Throughout Part II we use the symbol ‘A to denote 
the reflectance in air of a medium whose matrix may 
have a refractive index exceeding unity. The reflec- 
tance of a medium having the same values of m» and # 
but a matrix of refractive index unity is denoted, as 
before, by R. With this notation, ‘ and R are identical 
if NV (matrix) 1. 


VII. Reflection and transmission at a surface of refrac- 
tive index V.— We shall first review relations between 
internal and external transmission at a surface of 
refractive index NV. 

Consider a constant temperature enclosure contain- 
ing regions of refractive index unity and JV. Distin- 
guish quantities in the two regions by (m) and (u) 
respectively, where m and yu are cosines of the angles 
between a light ray and the normal to the surface at 
the point of incidence measured in the media of refrac- 
tive index unity and WV respectively ; then the radia- 
tion transmitted from the first to the second region 
in solid angle dQ(m) is J(m) mt(m) d2(m), where J 
is the intensity and ¢ is the transmittance of the sur- 
face ; and this must equal the radiation transmitted 
in the reverse direction, so 


(Tih) 
But 


I(m) mt(m) dQ(m) = I(p) ut(u) dQ(z). 


iy 


m dQ(m) = N2 yp dQ(p), 


and in a constant temperature enclosure 


I(u) 
I — ees 
(m) = 58 
Hence it follows that 
(ied) (mm) = ta): 


Now suppose that uniformly diffused radiation of 
intensity / be incident on the surface from within a 
medium of refractive index NV, the transmitted flux is 


amit 


uw t(u) du. 


For rays incident at angles exceeding the critical 
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angle, ¢(y) is zero. At other angles (7.2) applies. Thus 


A wad ay ii 
x ff = —— m t(m) dm. 
NN a 0 
« 
r1 
But 2 | m t(m) dmis the transmittance of the surface 
e 0 


for diffused radiation incident from the region of 
refractive index unity, denoted by J. Then 


3) 
N? 


A 
wT = gah 


so that the transmittance of the surface for diffused 
radiation incident from within is J/V?. 

The value of J for an optically plane surface has 
been given by Watsu [6], from whose expression the 
following values have been computed : 


N J 
1.333 0.933594 
1.50 0.908222 
1.52 0.905293 


VIII. Uniformly diffusing surfaces. — Let flux 
= Fj be incident from air from any or all directions 
onto unit area of the surface of a medium whose matrix 
is of refractive index V. The surface is supposed to 
possess the hypothetical properties that the radia- 
tions reflected from and transmitted through it are 
both perfectly diffused, and the transmission and 
reflection properties of the surface are independent of 
the direction of incidence. Then the flux transmitted 
through the surface is x Fi J, the flux reflected at the 
surface is x Fi(1 — J). 

Let the outward flux just below the surface of the 
diffuser be =. Then of this x JJ/N? is transmitted, 
and x4 (1 J/N?) reflected back into the medium. 
Thus the total inward flux immediately within the 
surface is x Fi J + x J(1 — J/N?), and if Rp be the 
total reflectance of the medium if in contact with a 
transparent medium of the same refractive index, then 


nF 


Tits = rs te 4 > ? 
x Fid + x d(1— J/N?) 
whence 
a Fj J Ry 
(8.1) nrAR= 


1 — Ry(t— dN) 


Thus the reflectance i of the diffuser in air, which is 
the ratio of the reflected flux to the incident flux, is 


Sues 
N2441—R,(t—J/N%)} 


Expressions for Ry appropriate to various phase func- 
tions have been given in Part I. 
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IX. Specularly reflecting surfaces. — Let a colli- 
mated beam of radiation of flux « Fi(m) per unit area 
normal to the beam be incident on the specularly 
reflecting surface of a semi-infinite diffuser at an angle 
to the outward normal whose cosine is — m. Then the 
incident illumination is =F i(m)m. The direct flux 
immediately below the surface, measured on a plane 
normal to the beam, is 

nr Fi(m) ae i(m) 
Uo 


per unit area and proceeds in a direction up. 


As in the previous section, consider the outward 
radiation field just below the surface. The flux incident 
on unit area of the surface from the range of direc- 
tions du may be written « S(u) du, or 2 x [(p) uw du, 
where /(u) is the average value of the intensity taken 
around all orientations at fixed u. This flux is reflected 
internally with reflection coefficient r(u) and, together 
with that part of the externally incident flux which 
has passed through the surface, constitutes the inward 
radiation field just below the surface. 

The outward flux results from the diffuse reflection 
within the medium of the two components of the 
inward radiation field. Denote the directional reflec- 
tance of the diffuser, in contact with a transparent 
medium of the same refractive index as the matrix, 
by R(u, uo) where uw and vo are the directions of obser- 
vation and incidence respectively. A(u, wu») follows 
immediately from expressions for the reflected inten- 
sity /(u) in direction uw by the relation 


Tr (u) 
Fi( tt) Bo 


where zw Fi(uo)9 18 the incident flux in direction uy 
per unit area of surface. The outward flux is then 
given by 


mJ(u) du =2nt,udu + 2nrJ,(u) du 


(9.1) Rv, bo) = 


where /, and /, are the intensities due to the two 
reflected components, 1.¢. 


nr J(u) du = 2x Fi(m) mt(m) R (u, uo) udu 


eit 


(9.2) 


+2nypdp | F(uo) (Uo) A(v, Ho) dio. 


w/ 0 


For some purposes if is more convenient to express 
(9.2) in terms of the intensity /(u) = (2, + Z,) rather 
than the flux : 


(9.3) L(u) = Film) m tm) R(y, bo) + 


il 
+2] He) rue) Ras to) to dt 


a 


The flux transmitted outwards through the surface 
is 


(9.4) 


a2 
nBp= 2m] Hu) wt(u) dy 


15 i, Ch 


oe) 


and the total reflectance ‘R is then obtained from the 
relation 


ip nF, > 
al = SSS ae ES, 
nm F'i(m) m 


(9.5) 


where the first term on the right hand side represents 
the diffuse reflectance and the second term the spe- 
cular reflectance. 


X. Scattering according to a phase _ function 
ol + 2 cos 6); Eddington’s approximation. — 
EppINGToN’s approximation has the merit that while 
it yields results of only limited accuracy, it does 
enable a general expression to be obtained for the 
reflectance, incorporating the albedo for single scat- 
tering and the parameter x in the phase function 
w(1 + x cos 6). 

By multiplying (9.3) by » and integrating with res- 
pect to yu, the equation is found 


| : 
(10.1) [ L(u) w du = > Film) m t(m) R(u9) + 
0) and 


e/ 


(4 


= | L (uo) 7(Ho) R (29) Uo Ato, 
0 
where 
(4 
Rly) = 2 fw lu, vo) de 
0 


a 


is the total reflectance for radiation incident in direc- 
tion uv, on the diffuser when in contact with a trans- 
parent medium of the same refractive index as the 
matrix. 

An approximate solution of (10.1) is obtained by 
assuming J(u) to be constant with pv. Then 


I = Fi(m) mt(m) R(wo) + 224, 


where 
1 
J — ( r(Uo) R( 9) Uo Ayo 
so that 
(10.2) aes Fi\(m) m We) Fi(uo) 


(a 


Since with the Eppincron approximation A(u,) is 
given by (4.1), equation (10.2) provides an analytic 
expression for / in terms of the phase function- 
mo(l + x cos 6). 

From (7.3), the flux transmitted outwards through 
the surface when /(y) is constant with p is 


xis 
N? 


80 that from (9.5) the total reflectance of the diffuser 
Is 
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. i(m) R(129) 
10.3 R = — + 
eee) N= 29) | 


XI. Solutions of higher accuracy. — The solution 
obtained above on the assumption of constant J(u) 
lies only part way towards the exact solution. To 
obtain improved accuracy, account must be taken of 
the variation of J(u.) with p, and this may be !done by 
replacing the integral in (9.3) by a quadrature for- 
mula. The most appropriate division of the interval 
of up from 0 to 1 offers considerable difficulty, for 
ideally a Gaussian division should be made. However, 
such a division depends on the weighting function 
(uo) R(u, vo) Yo, and while r(y») 1s small provided 
Uy > Ue, Where pv. is the cosine of the critical angle, 
it is just this part of the division which is important 
when (9.4) is to be evaluated subsequently. Further, 
R(u, vo) depends on both m, and yu, so that different 
divisions would be best for every mo and every p, 
which is clearly impracticable. Fortunately, the range 
of variation of R(u, wo) with vo is not very great when 
v) is not too small. The interval will therefore be divi- 
ded with the simple weighting function yw», in which 

af 
case | I (ue) r(eo) Alvi te) Uo deo im (9.3) may be 
al! @ 


replaced, for greatest convenience in subsequent 
computation, by » L(y.) R(w,, ) CF where according 


j 
to the method of Gaussian division, the weight a; 


attached to the division ; is 
4 4 = 
(11.4) a, = ———__ | ED ap | orto) de. 
Tl @—e)J 0 Lz; A 


Furthermore, the division points v.; are such that 


an n 


| f(u) udp = > C {(v,) 
0 


i=4 


(14.2) 


e 


where f(z) is any polynomial of degree 2n — 1 and the 
c, are appropriate numerical constants. Let 


Mu) = bo 4 


byt Oo? eee 


hen on substituting in (11.2) and equating coeffi- 
cients of bj, there follows the set of equations whose 
solutions give Bj: 


(41.3) y 
hea 


G. —_ 
rks hi 


r+ 2 


where in successive equations r takes the values from 
0 to 2n-1. Values of the division i for small values of 
n are 
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yj 
Ll ty = 0667 
De == (0) xO 
2) 
iy == UO) res 
i == OLR 
3 ee =O ROOT 


Ge == 096 


On making the above substitutions, (9.3) is replaced 
by the set of n simultaneous equations 


(11.4) Z(uy) = F(m) mt(m) R(u,,u,) + 


where in successive equations k takes the values from 
1 to n. The solution of the set of equations (11.4) 
yields /(u,) to the nth approximation. 

With the values of /(u,) solved, the flux passing 
out into the medium of incidence is given by (9.4), 
which is in turn given to the nth approximation by 


ii =e >; A 
j=1 


where the weight A; attached to the division t; is 


now 
BO Aa J, {IL = 
en 0 ties ea = 


The total reflectance follows from an application of 
(Jap). 

Since the effort involved increases greatly with 
the order of the approximation, it is desirable to 
know the accuracy achieved with various values of n. 
The first three approximations do show appreciable 
differences, as the following examples indicate. 


Uo H(U49) Apo. 


Isotropic diffuser, specularly reflecting surface, 


IS (Gonenoabe) == slays (Uy al 


any Isl approx. 2nd approx. 3rd approx. 
1.000 1.0151 1.0228 1.00019 
0.975 0.4671 0.4527 0.4685 


! 


The exact result for m — 1 is 1.00000, so that the 
third approximation is in error by only 2 parts in ten 
thousand. This does not automatically mean a higher 
approximation would necessarily give a better result, 
or that the third approximation is necessarily satis- 
factory for other values of mp, though it seems likely 
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that the accuracy will improve as the absorption 
increases. 

Two further checks have been obtained i carrying 
out the computation of the reflectance on the third 
approximation for (i) a polished medium of refractive 
index 1.333 and m, = 1, and (ii) the same medium 
in contact with a cover glass of refractive index 
1.523. The results are 1.00011 and 1.00012 as compa- 
red with the exact values 1.00000. 

Finally, the third approximation has been used for 
the case VN = 1; the results may be compared with 
the exact values given in Table I, and for a few values 
of wy are. 


Ty R (third approx.) R (exact) 


1.000 
0.975 
0.800 


0.99999 
0.64093 
0.28526 


1.00000 
0.64092 
0.28526 


The agreement is excellent, whereas the first ap- 
proximation gives an error of about 1 per cent in R. 

There are therefore good grounds for believing that, 
for N > 1 also, the third approximation will yield 
results of adequate accuracy for all photometric pur- 
poses. 


XII. Results. — Table 4 gives the total reflectance 
of a diffuser scattering isotropically, of refractive 
index 1.5, with a perfectly diffusing surface (equa- 
tion 8.2); the values of Ry have been taken from 
Table 1. 


TABLE 4 


Total reflectance for a diffuser with a perfectly diffusing 


surface. Isotropic scattering. N (matrix) = 1.5 

Ho Refleetance Bq Reflectance 
1.000 1.00000 0.80 0.2492 
0.999 0.8567 0.7 0.2028 
0.995 0.7235 0.6 0.1725 
0.990 0.6454 0.5 0.1506 
0.975 0.5269 0.4 0.1347 
0.950 0.4314 0.3 0.1203 
0.925 0.3756 0.2 0.1092 
0.900 0.3369 Q.1 0.0998 
0.85 0.2845 0 0.09178 


In Table 5 results are given for the total reflec- 
tances of diffusers scattering isotropically or according 
to a phase function m,(1 + cos 0), as derived on the 
third approximation. The surface is specularly reflec- 
ting and the matrix of refractive index 1.5. Values of 
diX/dNV, obtained by computing ‘f on the first approxi- 
mation for matrices of refractive index 1.50 and 1.52, 
are also included, these enabling the results to be 
extended over a useful range of refractive index. As 
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TABLE 9 


Total reflectances for isotropic scattering and for a phase 
function a(1 + cos 6), obtained on the third approxima- 
tion. Gradients dih/dN derived on the first approximation. 


N (matrix) = 1.5 ; specularly reflecting surface ; uy = 1. 
lsotropie To (4 + cos 0) 
ory — — === = = = = 
RK AR JAN R AR dN 
1.000 1.00000 0.0000 1.00000 0.0000 
0.999 0.8414 —0.162 0.8107 —().186 
0.995 0.6910 —().251 0.6392 —().263 
0.990 0.6020 —0.271 0.5418 —(.270 
0.975 0.4685 —0).260 0.4003 —(0),233 
0.950 0.3612 —().216 0.2936 —0.170 
0.925 0.3008 —().177 072355 —(0.123 
0.900 0.2600 —(0).145 0.1975 —().087 
0.8 0.1720 —0.058 0.1206 0.000 
0.7 0.1289 —0.006 O.0868 +().045 
0.6 0.1024 0.030 0.0683 + ().073 
0.5 0.0841 + 0.057 0.0571 +0.091 
0.4 0.0708 10.077 0.0500 +0,104 
0.3 0.0605 {0.094 0.0455 SIE {O)su tals! 
0.2 0.0523 10.107 0.0426 aE (Nal 2X0) 
O14 0.0456 + 0.119 0.0409 SJ Or125 
0 0.04000 +0.129 0.04000 +0,129 


indicated in § XI, there are grounds for believing that 
the errors in the reflectances derived on the third 
approximation are confined to the fourth decimal. 


TABLE 6 
Reflectances obtained with the Epp1iNGTon approximation 


for phase functions of the forms a(1 + x cos Q). 
N (matrix) = 1.5, specularly reflecting surface ; u) = 1. 


Isotropic (1 + cos 0) B_(1 + 0.5 cos O) 

Bo (a 

R Correction} {R Correction) 9 Corrected KR 

1.000 |1.0000 6.0000) 1.0000 0.0000) 1.0000 1.0000 
0.999 |0.8426 —0.0012|0.8137 —0.0030/0.8311 0.8290 
0.995 |0.6959 —0.0049|0.6483 —0.0091|0.6748 0.6678 
0.990 |0.6101 —0.0081|0.5532 —0.0114}|0.5850 0.5752 
0.975 |0.4808 —0.0123)0.4166 —0.0163]0.4520 0.4377 
0.950 |0.3780 —0.0168)|0.3130 0.0194] 0.3486 0.3305 
0.925 |0.3189 —0.0181/0.2561 —0.0206/ 0.2904 0.2710 
0.900 |0.2785 —0.0185/0.2185 —0.0210|0.2511 0.2313 
0.8 0.1890 —0.0170|0.1403 —0.0197| 0.1663 0.1479 
0.7 0.1431 —0.0142]0.1040 —0.0172]/0.1247 0.1090 
0.6 0.1141 —0.0117|0.0828 —0.0145]|0.0992 0.0861 
0.5 0.0934 —0.0093)| 0.0690 —0.0119}0.0817 0.0711 
0.4 0.0786 —0.0078|0.0594 —0.0094]0.0689 0.0603 
0.3 0.0655 —0.0050|0.0524 —0.0069/0.0591 0.0531 
0.2 0.0555 —0.0032|0.0472 —0.0046/0.0514 0.0475 
O11 0.0471 —0.0015|0.0431 —0.0022/0.0451 0.0433 
0 0.0400 0.0000) 0.0400 0.0000)}0.0400 0.0400 


Comparisons are made in Table 6 between the re- 
sults of Table 5 and those obtained for the same diffu- 
sers on the EppiINGToN approximation. The correc- 
tions required to results obtained on the EpprnGron 
approximation to yield those obtained on the more 
accurate third approximation are almost the same for 
x=0 and x =1, so that the simpler Eppincron 
approximation may be used for intermediate cases 
of the phase function m,(1 + x cos 6) and the appro- 
priate correction found by linear interpolation. Re- 
sults obtained in this way for = 0.5 (column 7) 
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should be in error by less than 1 in the third decimal 
place, and this has been confirmed by comparison 
with a few results obtained on the third approxima- 
tion. From the data of Tables 5 and 6, reflectances for 
other phase functions of the type mo(1 + x cos 8) 
may readily be obtained using three-point interpola- 
tion between the values of x. Further, it may be noted 
that to within + 0.001 the reflectance for «= 0.5 
divides the interval between reflectances for «= 0 
and 1, for the same mp, in the ratio 0.45 for both 
N =1 and 1.5, and the same division is therefore 
applicable for intermediate values of N. 


TABLE 7 


Total reflectances for isotropic scattering and for a phase 
function a (1 + cos 9), obtained with the third approxi- 
mation and with the EppInGron approximation. N (ma- 
trix) = 1.333, polished surface, 44 = 1. 


Third ° eye. Correction required 

Third | EDDINGTON to Epp1ineton 

approximation approximation approximation 

“9 Pied Le : = ane 

Iso- Hq (14 lso- =o (1+ ie ge Wp (1+ 

tropic | cos 9) trople | cos 9) Pome py cos 9) 
1.000 1.000 1.0000} 1.0000} 1.0000 0.0000 0.0000 
0.99 0.6520] 0.5934)! 0.6584! 0.6031 }—0.0064|—0.0097 
0.95 0.4073] 0.3324] 0.4229) 0.3512/—0.0156|—0.0188 
0.80 0.1898} 0.1271] 0.2082] 0.1486/—0.0184|/—0.0215 
0.60 0.1030} 0.0596} 0.1161|} 0.0768|—0.0131|—0.0172 
0) 0.0204} 0.0204} 0.0204) 0.0204 0.0000) 0.0000 


Table 7 contains reflectances computed on the third 
approximation for specularly reflecting media having 
N(matrix) = 1.333, and phase functions o)(1+2 cos 6), 
where « = 0 or 1; these are compared with results 
from the EppiInGron approximation, and it will be 
seen that the corrections required in the latter case are 
very nearly the same as those of Table 6. Corree- 
tions for intermediate values of N can readily be 
interpolated, and even without interpolation, cor- 


TABLE 8 


Total reflectances for isotropic scattering and for a phase 
function m (1 + cos Q), obtained on the third appro- 
ximation. N (matrix) = 1.333, in contact with a cover glass 
of refractive index 1.523 ; 5 = 1. 


Reflectance 
com - 
Isotropic To (1 + cos Q) 
1.000 1.00000 1.00000 
0.999 0.8684 0.8420 
0.995 0.7361 0.6886 
0.990 0.6541 0.5966 
0.975 0.5242 0.4559 
0.950 0.4155 0.3436 
0.925 0.3512 0.2799 
0.900 0.3067 0.2372 
0.8 0.2072 0.1476 
ORG 0.1565 0.1068 
0.6 0.1246 0.0838 
0.5 0.1021 0.0697 
0.4 0.0858 0.0607 
Ors 0.0729 0.0547 
0.2 0.0626 0.0509 
O.1 0.0541 0.0484 
0 0.04702 0.04702 
| 
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rections from Table 6 can be used to an accuracy nor- 
mally sufficient for any NV in the range 1.333 < N< 1.5. 
Table 8 contains reflectances obtained on the third 
approximation for a medium scattering isotropically 
or with a phase function m,(1 + cos 6), the matrix 
being of refractive index 1.333 and the medium in 
optical contact with a cover glass of refractive index 
1.523. The results in Table 8 apply specifically to 
diffusing media whose matrix is water, contained in a 
glass cell ; the refractive index of the latter is usually 
close to 1.523. Comparison with Table 7 shows that 
the cover glass has little effect for total reflectances 
above about 0.6. 


TABLE 9 


Directional reflectances for light incident along the normal. 
N (matrix) = 1.5 ; specularly reflecting surface. 
Isotropic scattering 


| 
Reflectance 
20 angle to normal 
ge NGS 30° 45° 60° Fe 
1.0000; 1.0164 1.016 1.014 1.0040 0.962 0.788 
0.999 | 0.8440 0.844 0.844 0.8381 0.806 0.661 
0.995 0.6814 0.682 0.683 0.6807 0.657 0.540 
0.990 0.5858 0.587 0.589 0.5870 0.568 0.468 
0.975 0.4437 0.445 0.447 0.4483 0.435 0.360 
| 0.950 0.3302 Orson 0.334 0.3360 0.328 0.272 
0.925 0.2670 0.267 0.270 0.2728 0.267 0,221 
| 0.900 0.2245 0.225 0.228 0.2301 0.225 0.188 
0.85 0.1691 0.169 0.172 0.1742 0.171 0.143 
0.80 0.1336 0.134 0.136 (Opa kas7/al 0.136 0.114 
O27, | 0.0896 0.090 0.091 0.0930 0.092 0.077 
0.6 0.0627 0.063 0.064 0.0653 0.064 0.054 
0.5 0.0443 0.044 0.045 0.0463 0.045 0.038 
0.4 0.0309 0.031 0.032 0.0323 0.032 0.027 
0:3 0.0205 0.021 0.021 0.0215 0.021 0.018 
0.2 0.0123 0.012 0.013 0.0129 0.013 0.011 
0.1 0.0056 0.006 0.006 0.0059 0.006 0.005 
TABLE 10 


Directional reflectances for light incident along the normal. 


N (matrix) = 1.5 ; specularly reflecting surface. 
Phase function : am(1 + cos 0) 
Reflectance 
Bo angle 10 normal 

0° ASS ove 45° 60° vow 
1.0000} 1.0164 1.016 1.014 1.0040 0.962 0.788 
0.999 | 0.8104 0.811 0.811 0.806 OM7 7G O63 
0.995 | 0.6248 0.626 0.628 0.627 0.606 0.499 
0.990 | 0.5201 0.521 05524 0.525 0.510 0.421 
OLG7omnOsS6Sou 0.570) Ors74> S037 7 0.368 0.306 
0:950 | 0.2559 0.257 0.261 0.265 0.262 0.219 
0.925 | 0.1950 0.196 0.200 0.204 0.202 0.170 
0.900 | 0.1555 0.157 0.160 0.164 0.164 0.139 
0.8 0.0765 0.078 0.081 0.084 0.086 0.073 
(O)s% 0.0426 0.043 0.046 0.049 0.051 0.044 
0.6 0.0245 0.025 0.027 0.029 0.032 0.028 
0.5 0.0139 0.014 0.016 0.018 0.019 0.018 
0.4 0.0075 0.008 0.009 0.010 0.012 0.012 
0.3 0.0036 0.004 0.005 0.006 0.007 0.006 
0.2 0.0014 0.002 0.002 0.0038 0.003 0.003 
0.1 0.0003 0.001 0.001 0.001 0.001 0.001 
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TABLE 11 


Directional reflectances for light incident along the normal. 
N (matrix) = 1.333, medium in contact with a cover glass, 
Ne eos 


Isotropic scattering 


Reflectance 
o angle to normal 
OF Kay? 30° 45° 60° 75° 
1.000 1.012 1.011 1.008 = 1.001 0.953 0.829 
0.999 | 0.865 0.866 0.866 0.863 0.825 0.720 
O5995) O72 20 0.722 0.724 0.724 0.695 0.609 
0.990 | 0.631 0.633 0.635 0.637 0.614 0.540 
0.975 0.491 0.493 0.496 0.501 0.485 0.428 
OIO0 Ors 77 O73) VOrS82) 09387 0.377 0.334 
0.925 | 0.308 9.310 0.314 0.320 0.312 0.278 
0.900 | 0.262 OFZ035) Be OSZO7— sOs2713 0.267 0.238 
0.8 0.160 0.161 0.164 0.168 0.166 0.149 
0.7 0.109 0.110 O411 0.114 0.113 0.102 
0.6 0.077 0.078 0.079 0.081 0.081 0.073 
0.5 0.054 0.054 0.056 0.057 0.057 0.052 
0.4 0,038 0.038 0.039 0.041 0.040 0.036 
0.3 0.025 0.026 0.027 0.027 0.027 0.025 
0.2 0.015 0.015 0.016 0.016 0.016 0.015 
0.1 0.008 0.008 0.008 0.008 0.008 0.007 
TABLE 12 


Directional reflectances for light incident along the normal. 
N (matrix) = 1.333, medium in contact with a cover glass, 
INGE O23: 


Phase function : a(1 -+ cos @) 


Reflectance 
% angle to normal 
0° 159 30° 45° 60° We 
1.000 1.012 1.011 1.008 1.001 0.953 0.829 
0.999 0.836 0.837 0.837 0.835 0.799 0.698 
0.995 | 0.668 0.669 0.672 0.674 0.649 0.570 
0.990 | 0.567 0.568 0.572 0.577 0.558 0.492 
0.975 0.415 0.417 0.423 0.430 0.419 0.373 
0.950 0.295 0.297 0.304 0.311 0.307 0.275 
0.925 | 0.229 0.231 0.237 0.245 0.243 0.219 
0.900 | 0.183 0.186 0.191 0.200 0.201 0.182 
0.8 1 0.092 0.094 0.098 0.105 0.108 0.101 
0.7 0.052 0.053 0.057 0.062 0.067 0.063 
0.6 0.030 0.031 0.034 0.039 0.043 0.041 
0.5 0,017 0.018 0.020 0.024 0.027 0.026 
0.4 0.009 0.010 0.011 0.015 0.017 0.017 
(033 0.004 0.005 0.006 0.008 0.010 0.011 


In all the above results the total reflectances include 
the specular component, whose magnitude is the 
reflectance for mw) — 0. Diffuse reflectances are rea- 
dily obtained by subtraction. Further, the values for 
m) — 1 are the exact values, rather than the slightly 
different values obtained on the third approxima- 
tion. 

Directional reflectances for normally incident radia- 
tion, as derived from the third approximation, for 
various media with specularly reflecting surfaces are 
given in Tables 9-12. In a few cases results are quoted 
to four decimal places, and these have been computed 
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using (11.4). Results given to three decimals have 
been obtained using graphical interpolation between 
the Z(u;) ; comparison with the computed values has 
shown that the results from graphical interpolation 
should be reliable to one unit in the third decimal 
place. In obtaining reflectances, account has been 
taken of the transmittance of the surface both for 
the incident and the reflected radiation. For the direc- 
tional reflectance along the normal, the specular com- 
ponent has not been included. A few selected compu- 
tations for these cases have indicated that to within 
+ 0.001 the directional reflectance for « = 0.5 divides 
the interval between those for « = 0 and 1, for the 
same mo, in the ratio 0.45. 
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Recent advances in Oplics, par E. H. Linroor, 
un volume relié toile de x-286 pages (23,5 x 15 cm), 
illustré de 113 figures et 8 planches hors-texte, 
Oxford University Press, 1955 (50 s). 


Cet ouvrage traile essentiellement des progres effec- 
tués récemment dans la théorie de la formation des 
images opltiques et dans la conception des instruments 
astronomiques. Il comprend quatre parties. 

Dans la premiére « The optical image », Vauteur étudie 
dabord l’aspect géométrique de la formation des 
images. 

Apres avoir donné (dans un formalisme nouveau ou 
Von utilise des quantités complexes) une classification 
des termes d’aberrations, l’auteur s’attache a recher- 
cher une méthode de réduction optimum des aberra- 
tions, basée sur la réduction du rayon de giration de 
Vimage géométrique. 

On passe ensuite al’étude de la tache de diffraction 
en présence d’aberrations et une documentation assez 
abondante est fournie sur les résultats obtenus par 
divers auteurs (cas d’aberration sphérique, de coma, 
WVastigmatisme ou de superposition). On aborde ensuite 
le probleme de la cohérence partielle par une méthode 
tres voisine de celle de Hopxrns et les relations de base 
sont établies. Le contraste de phase est enfin considéré 
comme un cas particulier. 

Dans la seconde partie « The Foucault test », Vauteur 
présente une théorie détaillée de la méthode de 
FoucauLtT: basée sur des considérations géométriques, 
cette méthode met manifestement en jeu des phénomeénes 


de diffraction. L’auteur justifie en optique ondulatoire 
les indications de Voptique géométrique, en mettant 
néanmoins en évidence l’influence de la diffraction (par 
exemple, la présence de franges brillantes au voisinage 
de discontinuités et aussi en discutant les erreurs 
systématiques qui peuvent résulter de l’interprétation 
purement géométrique de figures de Foucault en 
présence de faibles aberrations zonales). 

La troisieme partie « The Schmidt camera» est 
consacrée a la chambre de ScumipT et a ses dérivées. 
On étudie par la méthode des « spots diagrams » les 
aberrations géométriques de Vobjectif original et l’on 
propose diverses solutions pour équilibrer les aberra- 
tions sur la surface image sphérique et pour obtenir 
un champ plan. 

Enfin, dans la quatri¢me partie « Plate diagramm 
analysis and its applications », Vauteur étudie la 
méthode d’analyse des aberrations du troisiéme ordre 
(due a Burcu) ot les aberrations de chaque dioptre sont 
représentées par une lame déformée fictive située dans 
le plan du centre; il Vapplique ensuite au cas de 
Vobjectif de CAssSEGRAIN et de ScumiIpT-CASSEGRAIN 
(proposé par BAKER et par l’auteur). 

L’ensemble apporte une riche documentation et 
auteur nous fait pleinement profiter de son expérience 
sur la résolution mathématique de nombreux problémes 
fondamentaux d’optique instrumentale. L’ouvrage est 
abondamment illustré (en particulier de documents 
relatifs aux taches de diffraction en présence d’aberra- 
tions) ce qui en facilite beaucoup la 
(A. MARECHAL). 
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Radiative transfer in two dimensions 


J. T. JEFFERIES 


Division of Physics, National Standards Laboratory, 
Commonwealth Scientific and Industrial Research Organization, Sydney. 


SUMMARY. — A first approxunation is obtained for the distribution of inlensity in a semi infinite scallering medium illuminated 
by normally incident radiation whose flux is an arbitrary function of one Cartesian coordinate in the boundary plane. 


SOMMAIRE. — Une premiére approximation est obtenue pour la distribution de Vintensité dans un milieu diffusant éclairé sous V inci- 
dence normale par un rayonnement dont le flux est une fonction arbitraire de Vune des coordonnées cartésiennes dans le plan 


frontiére. 


ZUSAMMENFASSUNG. 


Es wird eine erste Ndherung fiir die Intensitdlsverteilung in dem Halbraum eines streuenden Mediums unter 


folgenden Annahmen aufgestellt : Die Strahlung falll senkrecht ein, wobei der Lichtstrom eine willkiirliche Funktion einer Koordi- 


nate der Grenzfldche ist. 


1. Introduction. — The diffusion of radiation 
through a scattering medium is described by the equa- 
tion of transfer whose solution gives the specific inten- 
sity J in the medium as a function of position and 
direction. Attention, so far, has been mainly confined 
to plane parallel media for which the sources of radia- 
tion and the physical properties are uniform over any 
plane parallel to the bounding surface. In these cases, 
clearly, the space variation of / can be described in 
terms of a single coordinate, namely the depth in the 
medium. 

Among the problems of this one dimensional type 
which have been solved, is that of the intensity dis- 
tribution in a plane semi-infinite scattering medium 
uniformly illuminated by normally incident radiation. 
The main purpose of this paper is to derive a first 
approximation to the solution of the transfer equation 
in the corresponding two dimensional case, 1. e. when 
the normally incident flux has arbitrary variation 
along one Cartesian axis in the boundary plane and 
is constant along the other. 

Two dimensional transfer problems occur naturally 
in many connexions, an example is in the discussion 
of the resolving power of a photographic emulsion 
where sideways scattering of radiation is of prime 
importance. 


2. The equation of transfer. — For an isotropically 
scattering medium the equation of transfer may be 
written (CHANDRASEKHAR [1]) 


(1) a a pe | Lao 
Ne 


x ds 4. 


where J is the specific intensity in the direction S, mo 
is the albedo for single scattering and x 1s the atte- 
nuation or extinction coefficient. In a Cartesian coor- 
dinate system (1) becomes 
(@) = [12 m2 on | =F [re 
M1. OL oy one An 


oe 


where /, m and n are direction cosines of J and the 
integral is taken over a solid angle 4 x. Equation (2) 
is the basic equation of the problem whose solution 
gives J as a function of position and direction. 

A number of methods for the approximate solution 
of (2) in the one dimensional case have been given 
(see e. g. KoURGANOFF [2]). In the method of discrete 
ordinates, which will be used here, the radiation field 
is divided into streams in certain specified directions 
and the intensity in each of these directions determin- 
ed. The points of division of the radiation field 
having been made, the integral in (2) is replaced by a 
suitably weighted sum over the intensities in these 
directions. As the number of points of division is made 
greater and greater, the calculated intensity in any 
direction may be expected to tend to a solution of (2). 
If the radiation field is axially symmetric, the inte- 
gration over the azimuth can be carried out at once 
and only an approximation to the integral over the 
polar angle is required. In the problems discussed 
here, J depends on both polar and azimuthal angles 
0 and ¢» and an approximation to the double integral 
over these angles is required. 

In two-dimensional problems, suppose J to be inde- 
pendent of the y coordinate and divide the radiation 
field into beams with azimuths — measured from the 
x axis — of + 7/4 and + 3 r/4. The approximation to 
the integral over © then becomes 


) 


270 ‘e 
| Pie eS 
Jo = 


For the division of 0, the angle which the beam 
makes with the z axis, we find it convenient to make 
a division into the directions uj = cos0; =+ 1,+¢ 
where ¢ is small. It will be noticed that the beams in 
the directions » = + 1 have an unspecified azimuth ; 
however, each may be regarded as the resultant of 
four degenerate beams with azimuths + 2/4, + 3 2/4. 
The purpose of the division into + ¢ instead of simply 
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Fia. 1. — Directions of the beams into which the radiation field is divided 


taking »=0 is to enable us to formulate correct 
boundary conditions before taking the limit ¢—> 0. 
Attempts to treat a division with any other values of 
uj have led to considerable mathematical complexities. 

For any particular azimuth 9, it then follows that 
with the above division for p, 

eo! a 

Ji, Me due SMe wo 

= ee 


where 6; is 1/3 for uj =-+ 1 and 2/3 for uj; = + «. 


The approximation to the integral over the solid angle 
may thus be written, 


O27 (eA a | 
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the summations being over the four directions of 
gj. From the symmetry of the problem we have 
1(9i, wi) = 1(— ¢i, wi) so that the integral in (2) may 
be written 


i=4 
The directions of the beams are shown in figure 1 where 
the beams of equal intensities in the directions + 9; 
are indicated by the same number. 


3. Collimated normally incident radiation. — Con- 
sider a plane parallel isotropically scattering medium 


illuminated by normally incident radiation of flux 
F(x) per unit area. The transfer equation for the diffuse 
radiation in the medium may he written, according to 
the above approximation, 


6) 2 


where uw = xz, ¢ = xz and the last term arises from 
scattering of the incident flux which has penetrated 
to a depth z. We try to solve the set (5) by a Fourier 
transform method. 

The Fourrer transform with respect to u of equa- 
tions (5) may be written 


Ay fi ; lo} w ; 
re = Ls) SS eee 
where 

1 Paes 
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\/2 Tw a 12) 
and 
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y(8) = Gin | eS” Fu) du. 


Substituting for n; and J;, six first order simulta- 
neous differential equations are obtained whose solu- 
tions may be found in the usual manner. For the 
homogeneous part of (6), 


1) = fe’ + gee + hy eter 


where the f;, g; and /; are functions of the parameter & 
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alone and the /’s are obtained as solutions of the cubic 
equation in k?, 


(8) “|e Eee aes 
Y 1—h? {—(y+ ek)? 1—(y—ek)?_ 


where 7 = i él ‘Qand y = w,/3. This equation has 
roots, to the first order in ¢«, 
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The negative roots for & are rejected as incompa- 
tible with the requirement that Ij; ->0 as ¢—>o. 

The particular integral is found by assuming a solu- 
tion of the form J; = p; ee e “&(2) whence it follows, 

+ 7 

on substitution in (6) that pj = 0,11, py = —2/¥. 

The integration constants f;, g; and h; are found 
ome ce coagtions. 2, — I, == I, = 0) at v= 0: 
These give 


| fa + 81 + Ay = (3/2 7) (8) 
ea, 

Fees 

( fs+ go + hs =0 


which are equivalent to the equations 
0-8 
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(10) 


J (u,o) 
0-4 
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where s(f) = < Sf; and s(g) and s(h) are similarly 
_) 
defined. 


Now letting <« > 0, two of the k’s, say k, and ks, 
tend to infinity in such a way that ¢k, and ek; remain 
finite, being given by (9). It follows then, from (11) 
that 


iS 
s(f) = 7 ae (1—k;,) v(8) 
3 (1—k,) | 
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where ¢j = « kj. 
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Fic. 2. — The variation of intensity at the surface for ™) = 1 and 0.75 with an incident flux 
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Fic, 3, — The variation of intensity at the surface for ™) = {£ and 0.75 with an incident flux 
. 2<1 
LEXQHD) = \ 1, | a | \ . 


Rather than deal with the separate intensities Ji, 
we shall consider the mean intensity defined by the 
equation 


At 9 =0, the Fourrer transform J of J may be 
written, 


(13) JE,0)= =e b= — [s(f)+s(g)+5(h)] —U2) 


which becomes, using (12), 


On inversion, (14) then gives 


1 pte » 
(15)  J(w,0) = ae | JE,Oye— iE uge 
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which may be evaluated numerically or by expan- 
ding (14) in a power series in &? and integrating term 
by term. In the latter case, 

12% I /112 
“Cora 
(2 7) 32a, n=O \n/ 
v(zZ)e—! Sue 
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(16) J(w,0) = — F(u) + 


at e 2 . : 
where ( aa is the binomial coefficient. 
n/ 

As illustrations of the solution (15), cases have 
been considered corresponding to the input functions 


(i) F(u) 


—e- July 5 
(ii) F(u) = | 


ib w\VY2<1 
0,|u yy Qe 


and results for m») = 1 and 0.75 are shown in figures 2 
and 3. As the series (16) converges rather slowly, the 
results were obtained by numerical integration. The 
discontinuities in figure 3 are results of the approxi- 
mation. 


The intensity inside the medium follows from the 
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condition that the solution for ¢ > 0 must fit smoothly 
onto that for ¢ = 0. It follows therefore that 


(17) w | Gara D arta (62 (0) hye 


and J (i, 7) is then found by inversion of (17). 


4. The reflectance. — An indication of the accuracy 
of the above may be obtained on comparing the 
reflectance found on this approximation with the 
exact value obtained by Giovanetut [3] from CHan- 
DRASEKHAR’S results. 

The outward flux of radiation per unit area Fy is 
given, in the limit «> 0, by 


27 
Q- 


it — Lom 


This may be written, in terms of the Fourrer trans- 
forms, 


2% 
FF; = a ile 
2x s(f) 
3 1 —+ ky 
which, from (12) becomes 
1—k 
18 i eee tay 
( ) 0 en ky } ( 4) 


The total energy passing outwards per second across 
an infinite strip of the surface of unit width and with 
its axis in the z direction is given by 


(19) ine F,du = 
ee te). te (Lh, 
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Similarly the total energy per second passing into the 
medium in a similar strip is 
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Using Fourter’s integral theorem, i.e. that 


+20 
ff f(q) 7% dy = 2x f(0), 
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the integrals (19) and (20) may be evaluated and the 
reflectance R obtained as the ratio 
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Values of R computed from (21) are compared with 
exact values in Table 1. The agreement, while not 
perfect, indicates that the accuracy of this approxti- 
mation is probably adequate for many purposes. 


TABLE 1 


Reflectance as a function of a 


Ho Equ (21) Exact value * 
0 0.0000 0.0000 
0.2 0.0200 0.0352 
0.4 0.0501 0.0834 
0.6 0.1010 0.1554 
0.8 0.2087 0.2853 
0.9 0.3333 0.4150 
0.95 0.4606 0.5356 
0.99 0.7072 0.7528 
0.999 0.8962 0.9128 
1.000 1.0000 1.0000 

* As given by GIOVANELLI (1955). 

5, Discussion. — The case of normally incident 


radiation treated here in a first approximation is 
probably the simplest which could be’ chosen. For 
angles of incidence other than normal and for finer 
divisions of the angular coordinate u, the above me- 
thod yields very complex integrals which suggest 
that a Fourrer transform method of solution would 
be of little use in extending the approximation. 
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Introduction. — The main reasons for re-examining 
the C.1. E. standard photometric and colorimetric 
data for the average eye were : 

(a) doubts about the correctness of the standard 
luminous efliciency function V5 in the ends of the 
spectrum, particularly the blue end, (b) small disere- 
pancies between the colour differences of certain tita- 
nium pigments, as observed and as computed from 
the C. 1. E. data, (c) the view of some workers in 
colour that colour measurements for a matching field 
larger than that used in determining the C. I. E. data 
might be more suitable as a basis for the standard 
system. To these reasons, may be added the remark 
that the method of determining the distribution 
coefficients, x, y, 2 by combining measurements of the 
unit coordinates for the spectral colours and a stan- 
dard white, with measurements of the V, function, 
is not an ideal procedure, when we consider that 
x“, Y, 2 are directly measurable without introducing 
any heterochromatic brightness-matching. This is not 
a new point of view. GuiLtp expressed a similar opi- 
nion in 1931 when the C. 1. E. colorimetric data were 
first agreed. He considered however that the method 
adopted was the best that could be done at the time. 

The N. P. L. investigation of colour-matching has 
aimed at a direct measurement of the distribution 
coefficients, entirely independent of heterochromatic 
brightness-matching. Both small and large matching 
fields are being studied. 

Knowing the distribution coeflicients, the unit coor- 
dinates are immediately calculable. But the V, fune- 
tion is not completely determined without some 
additional measurements involving heterochromatic 
photometry. Certain additional measurements of this 
kind have been made. It will be clearer if the question 
of the V, function is deferred until after the work on 
the distribution coefficients has been summarised. 

Because of delays for various reasons, it was not 
until August 1954 that the apparatus specially cons- 
tructed for this work was completed and calibrated. 
Since then a considerable number of measurements 
have been made but I must emphasise that the mean 
results about to be considered are preliminary ones 
obtained for a pilot group of 10 subjects. Several inte- 
resting if difficult issues are raised by the results for 


the pilot group, on which it has been necessary to 
make special measurements — measurements not yet 
completed. In the light of all the preliminary work 
and of the discussions at the present meeting, the 
final conditions for measurements on a much larger 
group of 50 or more subjects will be fixed. These it is 
hoped might be completed in about one year and 
could form a contribution to a revised set of standard 
data which the Commission might later adopt by a 
postal vote, perhaps after an intersessional meeting of 
colorimetric specialists. 


The N. P. L. trichromator. — A full description of 
the N. P. L. apparatus for this work was given re- 
cently (** Thomas Young Oration ”, 25th February, 
1955), but is not yet available in print. The main fea- 
tures are as follows. Three double ‘‘ monochromators ”’ 
of van Cirrert type with cancelling dispersion are 
mounted vertically one above the other. In the central 
spectrum of the middle tier, a single displaceable slit 
selects the monochromatic test colour which illumi- 
nates one half of the matching field. Three fixed slits 
in the upper tier select the monochromatic mixture 
primaries. These after recombination in the second 
stage of the monochromator provide a mixture which 
illuminates the other half of the matching field. The 
intensities of the primaries are independently variable 
by optical wedges inserted immediately after the slits 
in the central spectrum. The lower tier provides a 
similar mixture in variable amounts of three prima- 
ries — of the same wavelengths as the mixture pri- 
maries — which is added to the test colour and desa- 
turates the latter. The exit slits of all three tiers are 
imaged at a common point in space and the subject’s 
eye is brought into position so that the pupil centre 
coincides with this point and the whole of the light 
in the slit images — which all lie within a square area 
of about 2 mm side — is collected by the pupil. The 
subject sees the fields by the method of Maxwellian 
view. The photometric cubes are so made that the 
mixture (or comparison) beam from the top tier not 
only supplies the mixture half of the field but also 
provides a surround which extends the total area of 
illuminated field to a diameter of 14°. In order to esta- 
blish the psychological pattern of a bipartite field 
distinct from a surround, a circular, “ grey ’ demar- 
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cation line is introduced in the mixture beam. Cubes 
for two bipartite matching fields of diameters 2° and 
10° respectively with horizontal dividing line are used, 
two cubes for each field size being provided so that the 
test field can be arranged to occupy either the upper 
or the lower half-field. The controls of the mixture 
wedges (and of the desaturation wedges) are so arran- 
ged that the intensity of the mixture can be varied 
without altering the proportions of the components, 
or alternatively the components can be varied inde- 
pendently, or, finally, any two components can be 
varied keeping their relative proportions the same. 

The intensities in energy units of the test colour 
and of the primary stimuli were measured by means 
of a calibrated rubidium-on-silver photocell, perma- 
nently installed beyond the eye point of the apparatus 
so that frequent checks and corrections of the energy 
intensities emerging from the instrument could be 
made. 

The three primary stimuli used for most measure- 
ments were as follows : 


Wave-number Wavelength 
of centre A 


Overall 
band width 
15,418 cm-} 
18,997 
22,456 


FRE CEPrIMAT Vive se.) 
Green primary 
Binesprimary «=. - 


648.6 my, 555/23.1 em /my, 
526.4 512/14,2 
445.3 601/ 
The intensity level at which colour-matches were 
made can best be expressed by giving the luminance 
in photopic trolands of the comparison half of the 
field, the luminous efficiencies of the mixture prima- 
ries being assigned their C. 1. E. values. Evaluated 
in this way the levels used vary from one subject to 
another but the variation is considerable only in the 
blue end of the spectrum. Between 690 my and 490 muy. 
the troland value for a typical subject ranged between 
400 and 1500. It dropped to about 180 at 460 my and 
remained about that value down to 410 mu. Below 
410 my. and above 690 my the troland value dropped 
sharply to 10 at 392 mu and 50 at 730 my, for a typi- 
cal subject. The total band width of the test colour 
diminished steadily from 190 cm (10 my) at 730 my 
to 96 em-! (2.2 my) at 476 my, and then rose to about 
510 em! (8.0 my) in the extreme blue 400-392 mu. 
Some desaturation is unavoidable when matching 
spectrum colours but the minimum amount is used if a 
single primary is added to the test colour for desatu- 
ration and the mixture field contains only the other 
two primaries. Initially this procedure was adopted 
for all test colours except those at and in the imme- 
diate neighbourhood of the primaries, where small 
amounts of two desaturating primaries must be added 
and all three primaries must be represented on the 
mixture side. It is of course necessary to see that the 
subject can ‘“ bracket ’ the match i.e. he must be 
able to set each of the three variable primaries (mix- 
ture or desaturation) so that the amount of the pri- 
mary is (a) too great (b) too small, for a match to be 
possible by any adjustments of the other two variable 
primaries. Particular difficulties were experienced by 
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almost all subjects for test stimuli in the range 460 mu- 
420 mu (approximately, subjects vary somewhat in 
their ‘* worst ’’ wavelengths). The trouble can best be 
described as a tendency to blue/green degeneracy as 
a result of which the range on the green (or blue) pri- 
mary wedge over which ‘ acceptable ” matches can 
be made by adjustments of the blue (or green) and 
red primaries only, is unexpectedly large. The aecep- 
table wedge range for the green primary is particu- 
larly wide and this introduces an instrumental error 
associated with the bracketing method of setting a 
logarithmic control (Compare Trezona, 1953, 1954). 
The error from this cause was reduced by adding green 
primary to both sides, so that the matches of test 
colours below 460 my. were made in a less saturated 
colour field than would be possible theoretically with 
a blue primary of wavelength 445 mu. In fact, the 
saturation corresponded approximately to the opti- 
mum obtamable with a blue primary at 460 muy. 


Mean data for pilot group. — A series of colour- 
matches through the spectrum gives the quantities 
of the three instrumental primaries in mixtures which 
match measured energy intensities of test colours of 
different wavelength. For every wavelength one of the 
quantities is of course negative, i.c. the corresponding 
primary is mixed with the test colour. By assuming 
the validity of the additive law for colour-matching, 
the intensities a,, 6), ¢ of the primaries to match 
unit energy intensity of each spectral colour are obtai- 
ned. These data for each subject were immediately 
transformed to a set of reference primaries, which 
were spectral colours very close to the instrument 
primaries and located at wave numbers 15,500, 19,000, 
22,500 respectively (645.2, 526.3, 444.4 mu). As a 
result the colour-matching functions are obtained, the 
quantities x, y,, %, of the reference primaries to 
match unit energy intensity of the spectrum colours. 
The logarithms of the mean values of these quantities 
for the pilot group are given for both 2° and 10° mat- 
ching fields in Table 2 of the summary of results ap- 
pended. For comparison, the corresponding quanti- 
ties for the C. 1. KE. standard colorimetric observer 
and for the modified observer suggested by Jupp 
at the 1951 Meeting of the Commission, have been 
computed, and the values are given in Table 3 of 
the appendix. The transformations involved in such 
a comparison imply that the additive law in its broa- 
dest sense is valid for complete colour-matching. It 
will be recalled that BLorriau (1947) and TREZONA 
(1953, 1954) have obtained evidence that the additive 
law is not accurately obeyed, and this point will be 
considered later. Meanwhile additivity is taken for 
granted. 

There are very considerable individual variations 
in the colour-matching functions particularly in the 
extreme violet (400 mu and below) where subjects 
may differ by some two log units. It must be remem- 
bered that at least some of the factors which cause the 
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wide spread at the ends of the spectrum in the V5 
curves of a group of subjects will also operate for the 
colour-matching functions. 

Turning first to the mean 2° colour-matching func- 
tions, it may be said that for the ‘ red ”’ function the 
agreement with the C. I.E. values is fairly close 
except at wavelengths below 444 mu. Here the agree- 
ment is better with Jupp’s modified C. I. E. values 
although even then the new values diverge from 
Jupp’s in the violet. The new “ green ” colour-mat- 
ching function is materially higher than the C. I. E. 
and Jupp functions at wavelengths between the green 
and blue primaries, while for the new ‘* blue ” colour- 
matching function the situation is reversed. Both the 
‘oreen ”’ and ‘ blue ” funetions lie between the Jupp 
and C. 1. E. functions at wavelengths below the blue 
primary (444.4 my). 

To analyse further these differences, a derivation 
has been made of the unit coordinates in a WricutT 
system, that is, one in which the units of the red and 
green spectral primaries are adjusted to make equal 
the red and green components of a mixture matching 
a spectral-yellow (17 250 cm-1, 579.7 mu), while the 
unit of the blue spectral primary is adjusted to make 
equal the blue and green components of a mixture 
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the unit coordinates in the Wriacur system would be 
identical and the factor could be regarded as expres- 
sing a difference of transmission of a pigment layer in 
the eye, covering visual receptors of similar spectral 
sensitivities. In fact, the unit coordinates for the new 
2° measurements are quite close to the C. I. E. values 
for wavelengths 700 to 490 my, but below 490 my 
certain differences are in evidence. These have to be 
considered in relation to the individual variations. 
For the “ red ” unit coordinate, a plot of individual 
deviations from the mean suggests that the difference 
of the mean and the C. I. E. values may be significant 
in the range 430 to 470 my (fig. 1). If so the new data 
differ from the C. 1. E. in a way which cannot be 
attributed to differences of pigmentation of the groups 
of subjects concerned. 

The information given by the unit coordinates on 
the Wricut system can be supplemented by conside- 
ring the adjustments of the units of the primaries 
which have to be made to transform the data to this 
system. In other words we consider the relative values 
of the original colour-matching functions at the nor- 
malising wavelengths. These values are : 


Mean 2° C15. Mean 10° 


matching a spectral blue-green (20 500 cm-!, 487.8 my). i Zot ou cs BAS ERR eas 3.40 3.56 4.11 
(See. Table 1 of the appendix). If the new colour-mat- al 7h 
ching functions differed from the C.1. EK. only by a 99599 em 
factor dependent on wavelength in the same way for (487 8) hot ee oe 108 in 
all three, ‘‘ red ”’ ‘‘ green ” and ‘ blue ” functions, mu. 
+- : 
a (oaks ame i= aes al 
Shia INDIVIDUAL VARIATIONS FROM MEAN (PILOT GROUP): 
a RED UNIT COORDINATE 
AO EXTREME VALUES ie (W.DW. SYSTEM) 
+03 TWICE AVERAGE a 
(SIGNLESS) DEVIATION ! 
+-O2}+—© DEVIATION OF C.I.E. ¢ 
| 
+01 ° 5 
O 
©) 
—:-O2 
PRIMARIES 15,500 cm! 
= TOS 19,000 u 
22,500 » 
Tey NORMALISERS 17,500 cm7 
20,5005 
—-°O5 
25,000 22,500 20,000 17500 15000 12500 cm-! 
400 450 500 600 7OO mu 


WAVE NUMBER / WAVELENGTH 


Fie, 1. — Individual deviations from the mean for “ red” unit coordinate (2° field). 
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UNIT COORDINATES 
BASED ON 
EQUIENERGY WHITE 
w=} (Elle 
—- PILOT GROUP 
(2° FIELD) 


400 450 500 550 600 700 mp. 
WAVE NUMBER / WAVELENGTH 
Fie. 2. — Mean unit coordinates for 2° field 


compared with C. I. E. r, g, b values. 


At the yellow normalising wavelength the new 
2° data and the C. I. E. give red to green ratios (x,/y,) 
which are satisfactorily close. At the blue-green nor- 
maliser, the blue to green ratio (z,/y,) is considerably 
lower for the new 2° data than for the C. I. E., a diffe- 
rence which is formally consistent with a smaller ave- 
rage density of macular pigment for the present sub- 
jects (It may be observed that the ratio is still further 
reduced in the new 10° results). 

The modified C. I. E. observer proposed by Jupp 
in 1951 is sometimes referred to as an individual less 
pigmented than the original C. I. E. observer. But in 
so far as this description is correct, the pigmentation 
in question appears to be lens (or vitreous humour) 
and not macular pigmentation. Jupp’s modification 
has the effect in the main of raising the values of the 
colour-matching functions in the violet : it represents 
an observer with rather more macular pigment than 
the C. 1. E. observer. These conclusions follow from 
a consideration of the form of the spectral absorption 
curves of lens and macular pigments, and of the 
method used by Jupp to derive his modified func- 
tions. 

The apparent difference in macular pigmentation 
between the new 2° data and the C. I. E. cannot in the 
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main be attributed to the use of a possibly incorrect 
>, function in the derivation of the C. I. E. distribu- 
tion coeflicients. It represents a discrepancy with the 
Guitp-Wricht measurements themselves. This is 
shown by transforming the new data to the spectral 
primaries 700.0, 546.1, 435.8 my. and deriving new unit 
coordinates, after adjusting the units of the primaries 
so that equal quantities are required in a mixture to 
match the equienergy white (wavelength basis). The 
calculation involves the integration of the colour- 
matching functions over the equienergy spectrum. 
The unit coordinates so obtained are directly compa- 
rable with the r, g, b, unit coordinates (coefficients tri- 
chromatiques) tabulated in the C. I. E. 1931 Colori- 
metry Resolutions (fig. 2). The major differences 
occur in the blue green where, for example, new unit 
coordinates at 490 my. for the C. I. E. observer corres- 
pond approximately to those at 486.8 mu for the new 
29 data. The hue limen at 490 my. is of the order of 
1.0 mu, so that the difference between the C. 1. E. and 
the new 2° values is about three times what would be 
detectable by a particular subject. 

lor the present purpose a better way of assessing 
the discrepancy is to compare it with individual diffe- 
rences. In this comparison the individual data in 
terms of the original reference primaries — 15,500, 
19.000, 22,500 em-1 — were used. All the subjects had 
made matches with a special “* NPL bluish white ” 
at the same time as their spectrum colour matches. 
For each subject the quantities of the primaries to 
match this special white were taken as the units for 
the primaries, and the unit coordinates for a selection 
of wavelengths were then evaluated. For the C. I. E. 
observer, the derivation of the unit coordinates on 
this basis entailed a calculation from the relative spec- 
tral energy distribution of the special bluish white. The 
special bluish white was in fact obtained by removing 
the slit in the central spectrum of the middle tier of the 
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trichromator so that reconstituted white light of 
measured spectral energy distribution filled the test 
half of the matching field. A neutral and a blue filter 
of known spectral transmissions were inserted to pro- 
vide a white of suitable intensity and colour (fig. 3). 
In a rectangular chromaticity chart with the ‘* red ” 
and “ green ” unit coordinates as axes, the points cor- 
responding to the matches of a particular spectrum 
colour made by the 10 subjects of the pilot group form 
a ‘“‘ cluster ’’, near the centre of which the calculated 
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location of 579.7 my for C. 1. E. and for in- 
dividual observers of pilot group (2° field). 


point for the C. 1. E. observer should be found if the 
agreement is good. This is the case for longer wave- 
lengths as shown by the diagram for % = 579.7 
(17,250 cm-) (fig. 4). In the blue-green, however the 
C. 1. E. point lies just on or just beyond the ‘t edge ” 
of the cluster. This is illustrated by the plot for 
A = 487.8 (20,500 cm) (fig. 5). If the points corres- 
ponding to the two groups of subjects on which the 
C. I. E. observer is based (the Guitp and Wricut 
groups) could be shown in this diagram they would 
form clusters overlapping the present one to some 
extent. 

While it is possible that the average macular pig- 
mentation for the present group of subjects happens 
to be less than that of the Guinp-Wricur groups by 
an amount corresponding to the above difference, a 
‘“ sampling ” error of this magnitude seems impro- 
bable when account is taken of the individual spread. 
This fact and the discrepancy in the unit coordinates 
on the Wricur system, makes it necessary to look 
for other causes of the difference. But before doing so 
the large field results will be described. 

Two special questions arise in large field matching : 
possible intrusion of rod vision, and apparent colour 
inhomogeneities in a small area of the field round 
about the direction of vision. At the field intensities 
used for the present work it seems improbable that 
rod vision is modifying materially the large field colour- 
match settings except in the extreme red end and 
possibly in the violet. I think we can understand why 
that should be from our knowledge of the spectral 
sensitivity and threshold sensitivity of the rod as com- 
pared with the cone mechanisms. [ cannot enlarge on 
this point now and I shall merely indicate where in 
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the mean results the effect of the rod mechanism 1s, 
we believe, making itself felt. 

The colour inhomogeneities at the centre might be 
expected to make it difficult to do the actual colour- 
matching but they have in fact proved of little trouble‘ 
as subjeets soon learn to ignore the small area of the 
field where they occur. The inhomogeneities are no 
doubt the net result of local variations in retinal pig- 
mentation, in the spectral sensitivities of the under- 
lying end-organs and in the neural organization 
beyond them. They are visible only with certain 
rather special, highly saturated or dichroic stimuli and 
are, for example, hardly noticeable when matching 
the standard bluish white. 
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Fic. 5. — Chromaticity chart showing the 


location of 487.8 mu for C. I. E. and for indi- 
vidual observers of pilot group (2° field). 


The mean 10° colour-matching functions compared 
with the mean 2° functions show on the whole the 
kind of change to be expected from a reduction of the 
macular pigmentation. It appears however that the 
difference is not wholly explicable by assuming recep- 
tor systems with identical (or linearly related) spectral 
sensitivities under uniform pigment layers of diffe- 
rent densities. This is apparent by a comparison of 
the unit coordinates in the Wricut system. The nega- 
tive “‘ red” unit coordinate in the range 510 to 460 mu. 
is appreciably lower for the 10° data. The small drop 
in the “‘ red” and the rises in the “ blue ” and negative 
‘“ green ” coordinates for deep reds beyond 700 my. 
are probably the result of rod intrusion. From a plot 
of the ratios of the 10° to the 2° mean colour-matching 
functions it appears that for wavelengths below 
590 my. all three ratios have a maximum in the blue- 
green region recalling the absorption curve of macular 
pigment but the shapes of the curves are not identical 
as would be the case if the simple explanation men- 
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tioned earlier held good. There is also a difference in 
Shape for wavelengths greater than about 580 mu, 
which cannot be ascribed to a pigment and certainly 
not to the macular pigment that we know. At 700 muy, 
the ratio curves for the blue and green colour-mat- 
ching functions show sharp changes associated, it is 
believed, with rod vision. ; 

To return to the 2° results there are certain diffe- 
rences between the present observational conditions 
and those used by Wricut and by Guiry which may 
be the reason for the discrepancies with the C. I. E. 
already noted. The principal differences in observa- 
tional conditions with comments on their possible 
effects are given below. These comments are based on 
such auxiliary measurements as we have been able to 
make in the time available, but to deal satisfactorily 
with some of the possibilities more experimental work 
will be necessary. 


Intensity Level of Matches. — From measurements 
made at 700, 580 and 490 my it appears unlikely that 
differences in intensity level have played an important 
part. There is some evidence that in the spectral region 
420 to 470 mu, a higher precision of colour-matching 
might have been obtained in the new work by making 
the observation at a lower intensity level. 


Surround. — In the critical blue-green region cutt- 
ing out the surround has little if any effect on the 
match settings. In the deep blue and violet there is 
some indication of an effect but this requires confir- 
mation. 


Degree of saturation. — It has already been noted 
that the addition of more green primary than theore- 
tically necessary helped to stabilise the matches for 
wavelengths below 460 my. It does not appear howe- 
ver that, on the average, the resultant settings were 
much modified. Tests on the use of extra red desatu- 
ration in the blue-green region showed slight changes 
but not sufficient to account for the N. P. L./C. I. E. 
differences. 


Pupil conditions. — The effective entrance pupil 
had a slightly larger area than that used by Wricur 
but was considerably small than GuiLp’s (about one 
seventh the area). From what is known of the depen- 
dence of colour on pupil of entry, the resultant hue 
difference in the critical regions in the blue-green and 
yellow would hardly exceed 0.5 my, and this would be 
too small to remove the major discrepancy. It might 
be making a small contribution however. 


Different Choice of primaries. — The red primary 
is not materially different for all three investigations. 
GuiLp’s green primary had an effective wavelength 
at about 540 my compared with 530 (WricHT) and 
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526 (Stites and Burcu). Guitp and Wricut used a 
blue primary with effective wavelength at or near 
460 my whereas a blue primary at 445 mu was used 
in the present work. This last difference was thought 
to be the most likely source of the discrepancy, and 
for 3 subjects of the pilot group a series of 2° colour- 
matching measurements was made on the trichroma- 
tor with the blue primary shifted from 445 mp to 
460 mu. The observations were in the main confined to 
shorter wavelengths (526 my and below) where a diffe- 
rence would be most likely to appear. A similar series 
for the 10° field was made by another set of 3 subjects 
from the pilot group. The measurements were mainly 
confined to wavelengths below 526 mu. The results 
were transformed to the same reference primaries 
15,500, 19,000, 22,500 as before. The new colour- 
matching functions should be the same as the original 
ones for the subjects concerned if the additive law is 
valid. The mean curves — which with only 3 instead 
of 10 subjects are a good deal less smooth than the 
main data — show only a small and perhaps not signi- 
ficant difference. The slight effects are similar for the 
small and large fields. With the 460 blue primary the 
‘* blue ” colour-matching function in its principal 
spectral region — 526 to the violet end — is a little 
less sharp, while the ‘ green’ and ‘“ red ” functions in 
their main spectral regions — 460 my to 650 my and 
526 to the red end respectively — are very little alte- 
red. This has the effect that the z/y, ratio at 
20,500 .cm-? is brought a little nearer the C. I. E. 
value. The effect on the z/y, ratio at 20,500 cm 
(487.8 mu) of changing the blue instrument primary 
from 22,500 to 21,750 is shown by the following 
figures : 


z/y7, (reference primaries 15,500, 19,000, 22,500) 


Co Vestas, 2° 10° 
Instrument GUE 
Primary 22,500 21,750 22,500 Paes SAO) 
we noo 547 tail 1.08 


means for 
3 subjects 
(Ness 1, 7, 14) 


means for 
3 subjects 
(NIH il, 335 7) 


Although the ratio obtained with the 21,750 primary 
(2° field) moves somewhat nearer the C. I. E. value 
—and all 3 subjects show a change in the same direc- 
tion — the gap is not much reduced. Another compa- 
rison that can be made is of the unit coordinates in a 
Warrant system (fig. 6). For the 2° case there is per- 
haps a slight systematic increase in the amount of 
desaturating ‘ green ” primary for test colours of 
wavelength below 444 my. when changing the “ blue ” 
primary from 445 to 460 my. There is also some irre- 
gularity in the unit coordinates for the latter case at 
test colours round about 460 my. This may be connec- 
ted with the rather different matching methods which 
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* blue” primary from 445 mu. to 460 my. on mean 
unit coordinates for three observers (W. D. W. 
system); 2° field). 


must be used near the primary. For the 10° field, the 
unit coordinates obtained with the different instru- 
mental ‘‘ blue ” primaries show even less difference 
than for the 2° field. An entirely independent exami- 
nation of the effect on the unit coordinates of changing 
the instrumental ‘‘ blue ” primary from 445 to 460 mu 
has been made for his own eye by Professor W. D. 
Waricut using his well-known spectral colorimeter. 
He found that there was no appreciable difference in 
the unit coordinates for the two cases. 

To sum up, we have been unsuccessful so far in 
establishing that one particular difference in observa- 
tional conditions could account for the discrepancy 
between the Guitp-Wriacurt and the present 2° data. 
On the physical side, rechecks of the relative spectral 
sensitivity of our photocell, on which so much de- 
pends, have been very satisfactory and an independent 
determination of the colour of the bluish white light 
provided by the test tier of the apparatus, using a 
DonaLpson 6-stimuli colorimeter, shows no appre- 
ciable difference from the colour calculated from the 
measured spectral energy distribution. It seems that 
the discrepancy may be the accumulation of small 
effects connected with observational conditions, cou- 
pled with an unexpectedly large difference in ‘ effec- 
tive ” macular pigmentation of the groups of subjects 
concerned. Uncertainties on the latter account should 
be removed when the measurements are extended to 
a larger group of 50 or 60 subjects. 

But we shall have to decide before making these 
further measurements precisely what observational 
conditions we will employ. It will be most valuable 
to us and to others making similar measurements to 
hear the views of delegates to this meeting of the 
C. 1. E. on what the observational conditions should 
be if the results are to provide revised standard 
colour-matching functions. Perhaps the principal 
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question to be considered is the desirable size of the 
matching field. | may perhaps indicate some of the 
points of measurement which will arise. With a large 
matching field the nature of the surround is probably 
not critical. With a small matching field on the other 
hand, a dark surround, a surround of the same lumi- 
nance and colour as the matching field or, finally, a 
surround of the same luminance but white, may give 
some small systematic differences. The intensity level 
for large field measurements must be kept high to 
avoid rod intrusion, and there is in consequence less 
freedom in choosing “ pupil of entry ” conditions to 
approximate normal viewing. The choice of primaries, 
the degree of desaturation to be employed and whe- 
ther it should be a primary stimulus or perhaps a 
standard white, are questions which arise with any 
field size. All the points mentioned have to be consi- 
dered with the probable applications of the data in 
mind. A point of another kind is that a standardisa- 
tion on small field data would enable the valuable 
Guitp-Wricut measurements to be incorporated 
— with a little difficulty, it is true — in new standard 
colour-matching functions. 


Relative luminous efficiency function, Vx. — The 
V), function was originally defined to give information 
about the relative energy intensities of spectrum 
colours which produce the same sensation of bright- 
ness. In the derivation of the C. I. E. distribution 
coefficients or colour-matching functions from the 
Guitp-Wricut measurements of the unit coordinates 
of the spectral colours and a standard white, a Va 
function was used. The validity of this derivation 
depends on the V) function being a linear combina- 
tion of the colour-matching functions, and this would 
be true if the matching process used in deriving V) 
(a) obeyed the additive law, and (b) always assessed 
as equal fields in complete colour match. For this 
purpose it is not necessary for the Vq function to give 
correct information about the brighnesses of different 
colours seen in juxtaposition. If the colour-matching 
functions are directly determined, the Vx function 
ceases to be of primary interest in setting-up the 
standard colorimetric eye. However, it will probably 
be generally agreed that the V, function should con- 
tinue to be a linear combination of the colour-mat- 
ching functions, and once the latter have been fixed 
there will remain two independent constants, speci- 
fying the relative luminosities of the reference pri- 
maries, which can be assigned so as to provide a V 
function best suited to the requirements of photo- 
metry. 

In the N. P. L. work, the relative luminosities of 
the spectral reference primaries (15,500, 49,000, 
22,500 em-1) were determined by direct comparison 
with a ‘‘ white ” mixture matching the NPL bluish 
white. The mean relative luminosity factors so obtain- 
ed for the pilot group are compared with the C. I. E. 
values in the following Table : 
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rent. The resulting , functions are given in Table 4 

of the Appendix. Both synthetic 1) functions for the 400 450 500 mu 


2° case differ considerably from either the C. I. E. or 
the Jupp V , functions. 

A comparison of the synthetic V, functions, for 
the 2° and 10° fields based in each case on the colour- 
matching functions and luminosity factors for the 
appropriate field size is of interest. The effect of redu- 
ced macular pigment in the large field case is shown 
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function for pilot group in blue end of spectrum 
with two possible linear combinations of their 
colour-matching functions (2° field). 


by the higher values for the 10° field, in the spectral 
region centred on 22,000 em-? (455 muy). 

The N. P. L. work has also included a determina- 
tion of the Vj function for a 2° field in the blue end of 
the spectrum, by direct comparison with a mono- 
chromatic blue field of wavenumber 22,456 cm-! (the 
blue primary in fact). The mean results for a group 
of 28 subjects have been compared with the new 
Isoak-TEELE measurements and with the Jupp and 
C. 1. E. V, functions in the accompanying graph 
(fig. 7). Both the new direct determinations of V in 
the blue end of the spectrum confirm that the C. I. E. 
values are too low for wavelengths below 460 my. 
(21,750 cm-). The relative luminosities for wave- 
lengths in the blue end are in much better accord 
with Jupp’s modified V function. 

The 28 subjects used in the N. P. L. work on V) in 
the blue end included the pilot group of 10 subjects 
for whom the colour-matching functions had been 
measured. The question arises whether their mean 
direct measured V, curve in the blue end is a linear 
combination of their 2° colour-matching functions. 
Two linear combinations are compared with the 
directly measured Vq values in figure 8. The combi- 
nation using the subject’s own directly measured 
luminosity factors (against white) is in fairly good 
agreement with V) direct, but the values are relati- 
vely low in the blue green. The combination using 
Jupp’s luminosity factors deviates more widely and 
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in the opposite sense. There is little doubt that an 
intermediate choice of the luminosity factors would 
give a rather close agreement. But the possibility of 
fitting with a linear combination of the colour-mat- 
ching functions, a directly measured section of the 
Vx curve, confined to wavelengths below about 
480 my is not of particular significance. More interes- 
ting would be a similar test applied to the complete V) 
curve determined for the same subjects by some accep- 
ted procedure of heterochromatic photometry, flicker 
against white, for example. Sufficient measurements 
for such a comparison have not yet been made. 

A final point worthy of note is that the mean direct 
measurements of V, in the blue end for the pilot 
group of 10 subjects are in fairly close agreement 
with the means for the expanded group of 28 subjects 
at least if the two wavelengths below 400 my are 
excluded. It does not appear from this that the pilot 
group has exceptional pigmentation characteristics 
in the blue end. 
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The work described above has been carried out as 
part of the research programme of the National Phy- 
sical Laboratory, and this paper is published by per- 
mission of the Director of the Laboratory. 
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APPENDIX 


N. P. L. colour-matching investigation (1955) 


Mean results for pilot group of ten subjects 


W. S. Stites and J. M. Burcu 


Introduction. — The mean results are given for a 
pilot group of 10 subjects (ages, 20-53 ; mean age 31) 
who made measurements both for a 2° and 10° bipar- 
tite matching field with a surround extending to an 
outer diameter of 14°. The data are stated for refe- 
rence stimuli of wave-numbers (wavelengths) 
15,500 em (6,452 A) 19,000 (5,263) and 22,500 
(4,444), which are very close to the instrumental sti- 
muli (15,418 (6,486), 18,997 (5,264) and 22,456 (4,453), 
approximately). The unit coordinates are stated in a 
W. D. W. system, that is to say, the units of the 
colour-matching functions are adjusted so that the 
“red ” and ‘ green ” functions are equal at the nor- 
malising wave-number (wavelength) 17,250 (5,797) 
and the ‘*‘ green ” and ‘‘ blue ” functions are equal at 
the normalising wave-number 20,500 (4,878). The 
C. 1. KE. and Jupp’s modified C. 1. E. data — as given 
in the Secretariat Report of the Colorimetry Com- 
mittec(@. Li. Proc: Vol.) Balle (1951) shave 
been transformed to the same reference stimuli and 
normalising wave-numbers and are included for com- 
parison. Such a transformation assumes that the addi- 
tive law, in its broadest sense, is valid for complete 
colour-matching. 


On the view that the relative luminous efficiency 
function Vx should be a linear combination of the 
colour-matching functions, it is necessary only to 
determine the relative luminous efficiencies of the 
reference primaries, to derive V, from the colour- 
matching functions. The results of this derivation are 
given for both 2° and 10° fields using the subject’s 
own determinations of the relative luminous efficien- 
cies of the reference stimuli by direct comparison 
with white under the same field conditions. In addi- 
tion, for the 2° case, the subject’s mean colour-mat- 
ching functions have been combined using the relative 
luminous efficiencies of the reference stimuli taken 
from Jupp’s modified C. I. E. Vx functions (C. I. E. 
1951). The C. I.E. and Jupp V, functions are also 
included for comparison. 

In deriving the means for the group of 10 subjects 
for the following Tables, it is always the quantity 
itself — unit coordinate, colour-matching function, 
V;, function — which is averaged, not the logarithm, 
In the case of the colour-matching functions and the 
V; function, in default of making both the mean quan- 
tity and the logarithm of the mean, the latter has been 
chosen as perhaps more convenient in giving a quick 
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survey of the results. An alternative method of deri- 
ving mean unit coordinates is to take as basis the 
mean colour-matching functions and to proceed as 
though these represented the properties of a single 
‘mean ” subject. The values of the unit coordinates 
for this mean subject differ slightly from the means 
of the unit coordinates derived independently for the 
individual subjects, which are the means quoted in 
the Table. The difference rarely exceeds 0.001 except 
in the violet where it reaches .003 for one wave num- 
ber (10° field) and .004 for three wave numbers 
(2° field). 

The calculations of the unit coordinates and the 
log colour-matching functions were made to at least 
one more figure than the quoted values. The latter are 


INVESTIGATION OF 


COLOUR MATCHING Abhi 
rounded to the third decimal place. For the unit 
coordinates this has the result that the sum for a given 
colour may differ from unity by + 0.001 in some cases. 

The final Table gives the mean results of a direct 
determination of the V3 function in the blue end of 
the spectrum by the method of direct comparison 
with monochromatic blue of wave number 22,456 em-, 
A larger group of 28 subjects, including the pilot 
group, was used and the means for the whole group 
and for the pilot group only are shown separately. 
The Table also includes for comparison the mean 
results of a recent direct determination of the V 
function in the blue made by IsHak and TEELE using 
the method of flicker photometry with a blue compa- 
rison light. All these data are for 2° field. 


TABLE 1 


Unit Coordinates 


(W. D. W. System) 


Mean Results for Pilot Group and C. I. E. Values 

| 2° matehing field (Ga We Mh 10° matching field 

| 1/A x 15 500 19 000 22,500 15 500 19 000 22 500 15 500 19 000 22 500 
(cm-) (Aye | 

| 13 750 Dee it 12025 — 0,080 0.005 vey == (0)0ey7/ 0.001 0.999 — 0.023 0.024 
14 000 AS 1.029 SNUB 0.003 OBS == (0),0837/ 0.001 1.012 —0.028 0.016 
14 250 7018 1,030) —— 0.032 0.002 OSG — 020317 0.001 1.022 —=— 02031 0.009 
14 500 6 897 1.031 — 0.033 0.002 1.035 — 0.036 0.001 12029 ==02034 0.005 
14 750 6780 | 120297 == 03031 0.002 1.031 —0.032 0.001 1.023 — 0.028 0.005 
15 000 6 667 1.023 — 0.026 0.003 1025" == 02026 0.001 1.017 — 0.024 0.007 
15250 6 557 1.015 - 0.016 0.001 (ROG e e017 0.000 1.012 —0.016 0.004 

rela 500 6 452 1.000 0.000 000 1.000 000 000 1.000 000 000 

| 15 750 6 349 0.973 0.028 —0.001 0.972 0.029 —0.001 0.976 0.027 — 0.003 

| 16 000 6250 | 0.932 0.070 0.002 0.932 0.069 — 0.002 0.938 0.069  — 0.007 
16 250 G54 | 0.875 0.128 — 0.003 0.877 O12659—=01003 0.878 0.129 — 0.007 
16 500 6061 | 0.797 0.207 — 0.004 0.801 0.204 —0.005 0.804 0.203 — 0.007 
16 750 5 970 0.708 0.299 0.006 0.706 0.302 — 0.008 0.711 0.297 — 0.009 
17 000 5 882 0.604 0.405 0.009 0.601 OFF1O OF Oat 0.607 0.403  — 0.010 
17 250 5797 | 0.505; 0.505, — 0.011 0.507 0.507 —0.013 0.506 0.506 —0.011 
17 500 Te 0.401 0.613 —0.013 0.399 0.617 —0.016 0.412 0.601 — 0.013 
17 750 5 634 0.314 0.702 — 0.016 0.313 0.705 —0.018 0.317 0.698 — 0.016 

| 18 000 5 556 0.236 O78 ——OL017 0.236 0.783 — 0.019 0.239 0.777 —0.016 
18 250 5 479 0.166 0.851 —0.017 0.170 0.849 —0.019 0.170 0.847 — 0.017 
18 500 5 405 0.109 0.906 —0.016 0.109 0.907 —0.016 0.108 0.907 —0.015 
18 750 5 333 0.054 0.955 — 0.010 0.054 0.956- = 02011 0.052 0.958 — 0.011 
19 O00 5 263 OOO 1.001 —0.001 000 1.000 O00 0.002 0.999 000 
19 250 5 195 — (1.052 1.031 0.002 — 0.052 1.031 0.021 — (0.048 1.022 0.026 
19 500 By 1 Lass — 0.094 1.027 0.067 — 0.099 1.031 0.068 — 0.088 1.024 0.064 
19 750 5 063 — 0,135 0.987 0.147 — 0.135 0.986 0.149 — 0.122 0.984 0.138 
20 000 5 000 — 0.156 0.890 0.266 — ().1156 0.893 0.263 — 0.142 0.887 0.255 
20 250 4 938 — (0.161 0.745 0.416 === (aC) Os75i 0.412 — 0.138 0.735 0.403 
20 500 4 878 — ().145 ODT 0.573 =X 0a lays) (O) a7 7% 0.577 —— 0.124 0.562 0.562 
20 750 4819 — (0.112 0.405 0.707 ——() alo, 0.409 0.724 — 0.096 0.390 0.705 
21 000 4 762 — 0.083 0.270 0.812 — (ost 0.271 0.834 — 0.068 0,256 0.812 
21 250 4 706 — 0.056 0.168 0.888 — 0.075 0.172 0.903 — (0.045 0.155 0.890 
21 500 4 651 — ().041 Ont 0.929 =— (0) (aes 0.106 0.949 — 0.029 0.093 0.936 
21 750 4 598 — 0.025 0.062 0.963 ——=()037 0.066 0.970 — 0.018 0.052 0.966 
22 000 4 545 — 0.014 0.035 0.980 —= (OP? 0.038 0.984 — 0.010 0.026 0.984 
22 250 4 494 — 0.006 0.013 0.993 == 0.0010 0.017 0.993 — (0.005 0.012 0.993 
22 500 4 444 000: == 0.001 1.001 O00 000 1.000 O00 000 1.000 
22 750 4 396 0.003 — 0.007 1.003 0007  — 02012 1.001 0.004 — 0.008 1.005 
23 000 4 348 0.008 — 0.013 1.005 0.012 — 0.022 1.001 0.007 — 0.014 1.007 
23 250 4 301 0.012 — 0.018 1.006 OLOL6ES —=105029) 1.012 0.010 — 0.018 1.008 
23 500 4 255 0.015 — 0.021 1.006 0.019 — 0.033 1.014 0.014 — 0.021 1.007 
23 750 4211 0.020) == 0:022 1.002 0.022 — 0.037 1.015 0.018 Oe 0.024 1.007 
24 000 4 167 0.023 —0.022 0.998 0.023 — 0.039 1.016 0.020 — ().026 1.006 
24 250 4124 0.026 — 0.023 0.996 0.024 — 0.040 1.015 0.022 — 0.029 1.006 
24 500 4 082 0.029 —0.018 0.989 0.025 — 0.040 1.015 0.023 —— 0.027 1.003 
24 750 4 040 0.030 — 0.024 0.994 0.026 — 0.041 1.015 0.024 — 0.029 1.004 
25 000 4 000 0.031 — 0.024 0.993 0.027 — 0.041 1.014 0.025 —: 0.027 1.002 
25 250 3 960 0.030 — 0.027 0.997 0.027 — 0.041 1.014 0.025 — 0.029 1.004 
25 500 3 922 0.030 — 0.026 0.996 0.028 —0.041 1.014 0.025 — 0.030 1.005 
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TABLE. 2 


Log (Colour-matching Functions) 
Mean Results for Pilot Group 
(An “ N ” after the logarithm means that the function is negative) 


2° Matching Field 10° Matching Field 

1/2 a 15 500 19 000 22 500 15 500 19 000 22 500 
(cm-1) (A) 
13 750 7 OMB 3.591 5.529 (N) 6.675 3.593 5.369 (N) 5.076 
14 000 7143 2.001 5.972 (N) 6.889 AO 5.858 (N) 5305 
14 250 7 018 2.393 1.369 (N) Dp ili 2.416 4,303 (N) 5.455 
14 500 6 897 2.794 4.777 (N) 5.494 2.820 4.732 (N) 5.592 
14 750 6 780 1.144 3.096 (N) 5.808 teil 3.001 (N) 5.975 
15 000 6 667 1.466 3.346 (N) 4.305 1.489 3.258 (N) 4.421 
15 250 6 557 1.758 3.435 (N) iN PUTT iol 3.370 (N) 4.405 
15 500 6 452 0.000 6.844 5.576 (N) 0.000 5.835 (N) 5.650 
5) 750) 6 349 0.184 2.120 4.603 (N) 0.165 2.004 4.747 (N) 
16 000 6 250 0.322 2.676 3.038 (N) 0.306 Dyn 3.3171(N) 
16 250 6 154 0.410 1.050 3.348 (N) 0.402 2.967 3.426 (N) 
16 500 6 061 0.450 1.338 3.489 (N) 0.461 1.260 3.524 (N) 
16 750 5 970 0.450 1.549 3.759 (N) 0.471 1.488 3.671 (N) 
17 000 5 882 0.412 weil 3.915 (N) 0.451 1.668 3.766 (N) 
17 250 5 797 0.348 1.820 2.038 (N) 0.387 1.783 3.843 (N) 
17 500 5 714 0.255 1.911 2.125 (N) 0.320 1.880 3.942 (N) 
WA IXO) 5 634 0.135 1.955 2.182 (N) 0.204 1.942 2.015 (N) 
18 000 5 556 0.005 1.998 2.207 (N) 0.079 1.986 2.019 (N) 
18 250 5 479 1.839 0.020 2.194 (N) 1.914 0.007 2.013 (N) 
18 500 5 405 1.631 0.022 2.134 (N) 1.698 0.018 3.967 (N) 
18 750 5 Sis} 1,307 0.022 3.897 (N) 1.359 0.016 3.777 (N) 
19 000 5 263 4.906 (N) 0.000 4.911 (N) 3.782 0.000 4.074 (N) 
19 250 5 195 1.187 (N) 1.954 2.153 1.248 (N) 1.972 2.090 
19 500 5 128 1.370 (N) 1.880 2.571 1.468 (N) 1.928 2.434 
19 750 5 063 1.456 (N) 1.794 2.844 1.573 (N) 1.873 D730 
20 000 5 000 1.458 (N) 1.688 1.035 1.619 (N) 1.806 2.979 
20 250 4 938 1.449 (N) 1.585 1.202 1.618 (N) 12735 1.188 
20 500 4 878 1.415 (N) 1.483 1 SEP 1.620 (N) 1.667 1.378 
20 750 4819 1.375 (N) 1.402 52 1.576 (N) 1.578 1.546 
21 000 4 762 1.328 (N) 1.314 1.660 1.518 (N) 1.485 1.697 
21 250 4 706 1.238 (N) 1.186 1.782 1.429 (N) 1.361 1.831 
21 500 4 651 1.125 (N) 1.037 1.825 1.302 (N) 1.201 1.912 
21 750 4598 2.957 (N) 2.829 1.899 1.140 (N) 2.991 1.968 
22 000 4545 2.741 (N) 2.598 1.920 2.903 (N) 2.708 0.000 
22 250 4 494 2-384 (N) 2.196 1.970 2.540 (N) 2.348 1.994 
22 500 4 444 4,820 3.095 (N) 0.000 5.817 (N) 3.017 0.000 
22 750 4 396 2.209 3.958 (N) 1.981 2.437 2.181 (N) 1.971 
23 000 4 348 2.507 2.235 (N) 1.949 2.668 2.349 (N) 1.912 
23 250 4301 2.630 2.310 (N) 1.900 2.748 2.399 (N) 1.836 
23 500 4 255 2.702 2.338 (N) 1.842 2.789 2.376 (N) 1.754 
23 750 4211 2.705 2.252 (N) 1.745 2.807 2.346 (N) 1.662 
24 000 4 167 2.640 2.147 (N) 1.606 2.714 2.218 (N) TaD 
a ae ; RG ae 2.002 oY) tee one 2.125 (N) 1.364 
D4 750 4 040 9.206 7508 ae oe aor re ooo 
25 000 4 000 3.957 3.333 (N) 5 784 Ane re es ay 

0 3 3.33: De De ; (N) 2.691 

25 250 3.960 3.708 3.077 (N) 2.548 3.758 3.190 (N) 2.430 
25 500 31922 3.405 4.803 (N) 2.242 3.433 4.853 (N) 2.097 
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TABLE 3 


Log (Colour-matching Functions) 
. I. E. and Jupp’s Modified C. I. E. Values 


CA SNe? ae the logarithm means that the function is negative) 
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(Ce Jupp 
1/A r 15 500 19 000 22 500 15 500 19 000 
(cn (A) 

(approx. values : + .015. 
gaps values are not 
meovan (Cy Me TE), 

13 750 7 273 = = = == = 
14 000 7 143 2.082 4.079 (N) 6.614 — 
14 250 7 018 2.448 4.447 (N) 6.946 _ os 
14.500 6 897 2.814 4.803 (N) 5.287 ~~ 4.65 (N) 
14.750 6 780 1.174 3.114 (N) 5.591 -- 3.01 (N) 
15 000 6 667 1.484 3.335 (N) 5.816 = 3.16 (N) 
15 250 6 557 1.769 3.440(N) 5.929 = Sel7 CN) 
15 500 6 452 0.000 5.560 6.969 — = 
15 750 6 349 0.180 2.106 4.540 (N) — 2.24 
16 000 6 250 0.319 2.639 3.082 (N) — = 
16 250 6 154 | 0.410 1.016 3.464 (N) — == 
16 500 6 061 0.453 1.306 3.755 (N) — one 
16 750 5 970 0.455 W535 3.984 (N) a = 
17 000 5882 0.425 1.707 2.158 (N) = = 
17 250 5 797 0.369 1.817 2.267 (N) = = 
Wee 5714 0,285 1.922 2.373 (N) a = 
17 750 5 634 0.180 1.981 2.425 (N) = 4 
18 000 5 556 0.050 0.018 2.443 (N) = ae 
18 250 5 479 1.889 0.037 2.420 (N) = a 
18 500 5 405 1.674 0.041 2.319 (N) == = 
18 750 5333 | 1335 0.029 2.106 (N) = = 
19 000 5 263 0.000 6.952 (N) — = 
19 250 5 195 .20 204 (N) 1.948 2.287 = = 
19 500 5 128 1.398 (N) 1.864 D7 ila = a 
19 750 5 063 1.449 (N) 1.760 2.973 = ina 
20 000 5 000 1.437(N) ‘1.642 1.145 = 1.63 
20 250 4 938 7.419 (N) 1.531 1.303 = 1.52 
20 500 4 878 1.404 (N) 1.424 1457 = 141 
20 750 4819 1.388 (N) 1.326 1.607 — 1.27 
21 000 4 762 1.355 (N) Ae216 1.738 TSPX (ON9) 1.15 
21 250 4 706 1.278 (N) 1.087 1.840 1.23 (N) 2.98 
21 500 4 651 1.211 (N) 2.946 1.930 1.10 (N) 2.81 
21 750 £598 1.069 (N) 2.774 1.972 2.93 (N) 2.66 
22 000 4 545 2.854 (N) 2.543 1.992 2.69 (N) 2.42 
22 250 4494 2.505 (N) 2.188 1.998 2.33 (N) 2.02 
22 500 4444 — (5.591) 0.000 — = 
22 750 1 396 2 345 2.033 (N) 1.988 2.42 2.03 (N) 
23 000 1 348 2.560 2 252 (N) 1.955 2.61 2.29 (N) 
23 250 1301 2.619 2.311 (N) 1.893 2.67 2.36 (N) 
23 500 4 255 2 576 2.265(N) 1.780 2.73 2.35 (N) 
"23-750 4 211 2.458 2.136 (N) 1.608 Dan 2.32 (N) 
24 000 4 167 2.278 3.950 (N) 1.400 2.67 DO (UNI 
24 250 4124 2.096 3.756 (N) 1.196 Dewi 2.19 (N) 
24 500 4 082 3.869 3.517 (N) 2.950 2.45 2.07 (N) 
24 750 4 040 3 682 3.325 (N) 2.755 2.28 3.94 (N) 
25 000 4 000 O22 3.159 (N) 2.985 2.20 3.79 (N) 
5 250 3 960 3.309 4,940 (N) 2.364 2.02 3.63 (N) 
25 500 3 922 3.021 4,642 (N) 2.063 3.82 3.44 (N) 


22 500 
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TABLE 4 


Log (Relative Luminous Efficiency, VA) 


Means for Pilot Group and C. I. &. and Jupp’s Modified C. I. E. Values 


[Opt ACTA 


Linear combinations of colour-matching functions 
/ 5 2° (Subjeets’ 10° (Subjects’ 
V/A KN C. I. BE. Jupp (1951) own direct 2° (Luminosity own direct 
(em~™) (A) comparison factors from comparison 
luminosity Junpp’s 1951 luminosity 
factors) curve) factors) 
13 750 TE PATS 4.801 As for C. I. E. 4,929 4,730 4,959 
14 000 7 143 3.192 > 3.344 3¥, 11257) 3.381 
14 250 7 018 3.558 » Bets 3.529 3.780 
14 500 6 897 3.926 , 2.136 3.929 2.182 
14 750 6 780 2.288 » 2.487 2.281 2.534 
15 000 6 667 2.602 2.811 2.607 2.854 
15 250 6 557 2.894 1.105 2.907 1.152 
15 500 6 452 Helse ies? 1.161 1.366 
15 750 6 349 1.338 D 1.550 1.366 1.542 
16 000 6 250 1.506 » 1.707 1.538 1.701 
16 250 6 154 1.640 » 1.821 1.671 1.816 
16 500 6 061 ets » 1.906 ean 1.904 
16 750 5 970 1.826 » 1.957 1.853 1.961 
17 000 5 882 1.891 » 1.989 1.912 1.999 
17 250 57907 1.941 » 0.005 1.950 0.008 
17 500 5 714 1.974 » 0.015 1.982 0.021 
17 750 5 634 1.993 » 0.002 1.988 0.010 
18 000 5 556 0.000 , 0.000 0.000 0.000 
18 250 5 479 1.996 » 1.986 1.999 1.977 
18 500 5 405 1.981 » 1.959 1.984 1.951 
18 750 DESGS 1.954 » 1.936 1.970 1.918 
19 000 5 263 1.910 » 1.890 1.934 1.875 
19 250 5 195 1.845 » 1.826 1.876 1.821 
19 500 5 128 1.749 » 1.741 1.792 1-755 
19 750 5 063 1.636 » 1.645 1.696 1.677 
20 000 5 000 1.509 » 1.527 1.582 1.591 
20 250 4 938 1.391 > 1.414 1.473 sult 
20 500 4 878 Nay » 1.314 1.369 1.435 
20 750 4 819 ea » e239 1.296 1-355 
21 000 4 762 1.074 » e777 1223 1.292 
21 250 4 706 2.970 » 1.099 1.130 1.228 
21 500 4 651 PRA » 1.005 1.021 1.162 
21 750 4 598 Q774 » 2.921 2.910 1.093 
(Approx.) = 
22 000 4 545 2.672 DGD 2.840 2.808 1.037 
22 250 4 494 2.568 2.67 2.811 2.724 2.995 
22 500 4 444 2.463 2.63 2.789 2.654 2.963 
22 750 4 396 2.351 2.59 2.758 2.587 2.904 
23 000 4 348 2.220 Disa 2.708 2.489 2.841 
23 250 4 301 2.069 Deets 3.653 2.402 2.765 
23 500 4 255 3.886 2.36 2.593 2313 2.700 
23 750 AY SGI 3.661 Py ON7 2.528 2.246 2.621 
24 000 4 167 3.422 2.16 2.408 PE) 2.508 
24 250 4124 Rill 2.03 2.279 3.979 2.348 
24 500 4 082 4.966 3.92 2.123 3.859 2.200 
24 750 4 040 4.770 3.80 3.907 3.606 3.994 
eee He ee mse ee 
De 4. 3.47 3.415 3.126 3.500 
25 500 3 922 4,089 3.29 3.104 4.797 3.168 


vor. 2 n° 4, Dee. 1955 | 


REPORT ON N. P. L. INVESTIGATION OF COLOUR MATCHING 


TABLE 


Logi) (Relative Luminous Efficiency, VA) 
Direct Measurements for 2° Field 
(VA adjusted to have value unity at 21 750 em-1) 
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1/A dr N. P. L. Pilot Group N. P. L. Extended Group IsuaKk-TEELE 
(cm71) (A) 10 Subjects 28 Subjects 104 Subjects 
20 750 4819 0.357 Ons82 — 
20 790 4 810 = ae 0.509 
21 000 4 762 0.288 0.321 = 
Diez 4 720 — es 0.368 
21 250 4706 0.203 0.222 — 
21 500 4651 0.104 0.120 = 

| 21 750 4598 0.000 0.000 = 
21 786 4 590 =. a 1.995 
21 978 4 550 — — 1.918 
22 000 4 545 1.909 1.904 — 
22 250 4494 1.857 1.855 — 
22 371 4470 = — 1.796 
22 500 4 444 1.811 1813 — 
22 624 4 420 = 1.830 
22 750 4 396 1.767 1.765 
23 000 4 348 1.703 daz — 
23 250 4 301 1.631 BEG 3H — 
23 500 4 255 1.585 1.600 — 
23 750 4 211 1.547 1.549 -— 
23°753 4 210 = == 1.346 
24 000 4 167 1.397 1.428 Macs 
24 213 4 130 — — 2.982 
24 250 4124 1.265 1287 
24 500 4 082 By) 1.136 — 
24 750 4 040 2.893 2.910 — 
25 000 4 000 2.693 2.693 = 
25 126 3 980 = = 2.694 
25 250 3 960 2.407 2.465 
25 500 3 922 2.128 2.214 — 
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Lettre a l’editeur 


Compensateur biréfringent a grand champ 


M. Francon et B. SERGENT 


Institut d’Optique, Paris 


On emploie généralement le compensateur de Basr- 
Nev en lumiére paralléle pour la mesure d’une diffé- 
rence de marche produite par une lame cristalline. 
Si on Tutilise en lumiere méme peu convergente le 
champ est limité par les lignes isochromatiques qui 
sont des hyperboles : Vordre d’interférences varie 
dune fagon continue dans tout le champ et introduit 
des erreurs dans les observations. I] est possible d’évi- 
ter cet inconvénient au moyen des combinaisons sul- 
vantes utilisant des lames taillées parallélement a 
axe: 

1) Si on remplace une lame cristalline d’épaisseur e 


par deux lames identiques d’épaisseur — croisées en 


2 

interposant une lame demi-onde, les biréfringences 
des deux lames s’ajoutent tandis que leurs variations 
qui sont de signes contraires se compensent en grande 
partie. Les lignes isochromatiques deviennent des 
cercles et le champ obtenu est 26 fois plus grand que 
celui dune lame simple pour du quartz. 

2) On obtient un résultat analogue en remplagant 
une lame critalline par deux lames taillées Pune dans un 
cristal positif (quartz) autre taillée dans un cristal 
négatif (spath) le rapport des épaisseurs étant conve- 
nable. 

3) Un troisieme systeme consiste & superposer une 
lame taillée dans un premier cristal et 2 lames d’axes 
croisés taillées dans un deuxiéme cristal de signe con- 
traire a celui du premier. Le choix des épaisseurs étant 
détermine. 

Ces différents systemes ont été indiqués et utilisés 
par B. Lyor dans son filtre monochromatique polari- 
sant. Ils peuvent étre également employés pour réali- 
ser un compensateur du type BABiner a grand champ. 

En adoptant par exemple la premiére solution on a 
la figure 1. Les quatre lames A BC D sont quatre 
lames de quartz. Les lames A et D sont identiques et 
leur 6paisseur est égale a l’épaisseur moyenne des 
lames prismatiques B et C. Ces 2 lames sont égale- 
ment identiques et leur ensemble forme une lame A 
faces paralléles. Les axes des lames C et D sont tour- 
nés de 90° par rapport a ceux des lames A et B. Entre 
A et B et entre C et D sont intercalées 2 lames demi- 


ondes E et F dont les axes sont a 45° des axes des lames 
en quartz. Il peut étre intéressant dans certains cas 
de rendre les lames A et D également prismatiques 


j ; 
H ; 
H r] 
i ry 
§ Y 
4 j 
j y 
j 
i Y 
4 4 
3 F 

1guver 41 


pour conserver Pégalité des épaisseurs A et B, C et D. 
L’ensemble des lames restant assimilable a une lame 
a faces paralléles. 
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Le compensateur formé par ces lames montre des 
franges en lumiére blanche comme dans le compensa- 
teur de BaBineT ou le prisme de WoLLAsTon et locali- 
sées également sur la face inclinée commune aux lames 
B et C. Avec un BaBINET ou un Wo taston la diffé- 
rence de marche varie rapidement avec V’inclinaison 
des rayons. Sur la figure 2 le compensateur de Bant- 
NET Q (ou le prisme de WoLLAsTON) est traversé par 
un faisceau convergent de lumiére en formant l’image 
@une source étroite sur Q. Un objectif O placé der- 
riere Q montre dans son plan focal = les lignes isochro- 
matiques des lames taillées parallelement a Vaxe. Ce 
sont des hyperboles et le champ angulaire 0 corres- 
pondant a une différence de marche pratiquement 
constante c’est-a-dire a une teinte uniforme est petit. 
Si on remplace le Basinet ou le WoLLASTON par le 
nouveau compensateur on obtient un champ beaucoup 
plus grand et on peut utiliser un faisceau de lumiere 
fortement convergent. 
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de Smiru dans lequel on a remplacé les deux prismes 
de WoLLaston par deux compensateurs du nouveau 
type mais au lieu de disposer le compensateur dans le 
plan focal de lobjectif du microscope il est possible 
maintenant de le placer dans loculaire. L’image de la 
source lumineuse se forme en S au foyer du conden- 
seur C. La préparation est en P et son image donnée 
par lobjectif O, du microscope, en P’. Un objectif 
auxiliaire O, donne du plan focal S’, généralement non 
accessible, une image en S” sur laquelle se trouve le 
compensateur Q,. Le compensateur Q, possédant un 
grand champ il est possible de le placer en S” ot l’ou- 
verture des faisceaux est plus grande qu’en S’. Ceci 
ne serait pas possible avec un BABINET ou un WOLLAS- 
TON dont les champs sont tres petits. Quant 4 image P’ 
Pobjectif O, la transporte en P” ot elle est observée 
par loculaire O¢. Un deuxieéme compensateur Q, peut 
étre placé en S conjugué de S" a travers C, O, et Oy. 
la compensation est ainsi assurée et il est possible 


Le compensateur a grand champ peut étre employé 
en microscopie interférentielle ol grace a ses propriétés 
il est possible de le placer dans l’oculaire du micros- 
cope. 

Des montages interférentiels simples peuvent étre 
réalisés pour létude des défauts d’homogénéité ou 
des défauts de poli dans le domaine macroscopique. 
De larges surfaces peuvent étre étudiées. 

La figure 3 représente adaptation du compensa- 
teur a grand champ au microscope. C’est le dispositif 


dutiliser en S une source large. Si le faisceau n’est pas 
trop ouvert ce qui est en général le cas, un simple 
Wot taston suffira en S. Naturellement on place un 
polariseur avant Q, et un analyseur apres Q,. 

La propriété d’autocompensation du prisme de 
WoLLASTON traversé deux fois; a été montrée par 
Léon et Francois Lu NouvEL puis utilisée par NOMARSKI 
en miucroscopie. Cette propriété est conservée ici, 
L. et F. Lenouvet ont indiqué un montage par 
autocollimation trés lumineux quwil est également 
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possible utiliser. L’objet transparent macroscopique 
est en P contre un miroir sphérique M (fig. 4). Au 
moyen dun objectif O, et aprés réflexion sur un 
miroir G on forme en S’ une image de la source S. 
Liimage S’ est au voisinage du centre du miroir M 
qui en donne une nouvelle image en S". Le compen- 
sateur Q a grand champ est placé sur les images 
S’ et S". Un objectif O, donne de P une image en P’ 
qui est projetée sur un écran E ou observée a l’ocu- 
laire. Le polariseur P est placé avant Q et Pana- 
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lyseur A aprés Q sur le faisceau de retour. Grace 
au champ considérable du compensateur il est pos- 
sible d’utiliser un miroir M d’ouverture « impor- 
tante et d’observer des objets P trés étendus. 

Un montage pour l’observation des objets transpa- 
rents macroscopiques traversés une fois pourra étre 
réalisé sur un principe analogue a celui de la figure 3 
et avec deux objectifs seulement. 


Manuscerit recu le 20 seplembre 1955. 


N° d’éditeur 122, — Revue @ optique. 


— Le Direcleur ; Emile PERRIN. 


Commission paritaire des papiers de Presse no 29754 


1981. — Impr. JOUVE, 15, rue Racine sParisy=== (Oss 5 


| 
| 


DATE DUE | 
UEt Of PER 


NOV 10 1983 
ii: OV 0 1983 


Layee eal: papers ; 
TO PERIODICALD DEsx © 


DEMCO 38-297 


